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Abstract 

In this paper we venture a new look at the linear isotropic indeterminate couple stress model in the 
general framework of second gradient elasticity and we propose a new alternative formulation which obeys 
Cauchy-Boltzmann’s axiom of the symmetry of the force stress tensor. For this model we prove the existence 
of solutions for the equilibrium problem. Relations with other gradient elastic theories and the possibility 
to switch from a 4th order (gradient elastic) problem to a 2nd order micromorphic model are also discussed 
with a view of obtaining symmetric force-stress tensors. It is shown that the indeterminate couple stress 
model can be written entirely with symmetric force-stress and symmetric couple-stress. The difference of 
the alternative models rests in specifying traction boundary conditions of either rotational type or strain 
type. If rotational type boundary conditions are used in the partial integration, the classical anti-symmetric 
nonlocal force stress tensor formulation is obtained. Otherwise, the difference in both formulations is only a 
divergence-free second order stress field such that the field equations are the same, but the traction bound¬ 
ary conditions are different. For these results we employ a novel integrability condition, connecting the 
infinitesimal continuum rotation and the infinitesimal continuum strain. Moreover, we provide the com¬ 
plete, consistent traction boundary conditions for both models. 
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1 Introduction 


1.1 General viewpoint 

The Cosserat model is an extended continuum model which features independent degrees of rotation in addition 
to the standard translational degrees of particles, see [551IH31 ISIl ISHl IM] for a detailed exposition. The prize, 
which has to be paid for this extension are non-symmetric force stress tensors together with so-called couple 
stress tensors which then represent the response of the model due to spatially differing Cosserat rotations. The 
couple stress model is the Cosserat model [48] with restricted rotations, i.e. in which the Cosserat rotations 
coincide with the continuum rotations. As such it belongs also to a certain subclass of gradient elasticity 
model^, where the higher derivatives only act on the continuum rotations. This constitutes a big conceptual 
advantage since the interpretation of the Cosserat rotations as new physical degrees of freedom is in general a 
difficult task. Such a model is also called a model with “latent microstructure” [TTl [12] . 

Let F = RU he the polar decomposition of the deformation gradient F = Vcp into rotation R G SO(3) 
and positive definite symmetric right stretch tensor U = VF'^F, where C —!> characterizes the 

deformation of the material filling the domain 17 C We write R — polar(F). In a variational context, the 
energy density W to be minimized in the geometrically nonlinear constrained Cosserat model is given by 

W = W(U_^, polar(F)^VxPolar(A)), (1.1) 

strain curvature 

which reduced form follows from left-invariance of the Lagrangian W under superposed rotations. In this paper, 
our objectives are much more modest. We will only be concerned with the linearized variant of (inj, which 
can be written as 

W = W(symVu, V[axl(skewVM)]) = Wiin(symVu) -I- Wcurv(V[axl(skewVu)]), (1.2) 

infinitesimal infinitesimal 
strain curvature 

where u : 17 C R.^ —t is the displacement and 

V axl(skewVu) = 2curltt. (1-3) 

The energy density (HH) is the classical Lagrangian for the indeterminate couple stress formulation. As will be 
seen later, this formulation leads naturally to totally skew symmetric nonlocal force stress contributions. 

Toupin already remarked on an alternative representation of the energy (fTT]) fTU51 Section 6] which leads, 
in its linearized variant given by Mindlin [71 eq. (2.4)] to a dependence on 

W = ir(sym Vu, Curl (sym Vu)) = Wiin(symVu) -I- Wcurv(Curl (sym Vu)), (1.4) 

due to the equivalence 

V(axlskewVu) = -Vcurlu = (Curl (sym Vit))^ 

instead of (HU. The representation iVcurlu is directly derived from the original Cosserat model [UlinH]. Both 
authors, Toupin and Mindlin, noted that now, comparing and the force stress tensors and the couple 
stress tensors are changed while the balance of linear momentum equation remains unchanged such that these 
concepts are not uniquely defined (see also Truesdell and Toupin’s remark on null-tensors [1051 p. 547]. However, 
they apparently did not realize that it is possible to use this ambiguity to obtain completely symmetric force 
stress tensors also in the couple stress model which is otherwise the paragon for a model having non-symmetric 
force stress tensors. We also need to remark that in a purely mechanical context, the observation of size-effects 
does not necessitate to introduce skew-symmetric stress-tensors |5]. 

In this paper we do not discuss in detail the field of applications of such a special format of gradient elasticity 
model. Suffice it to say that much attention is directed to nano-scaled material in which size-effects may become 
important, which may make the present model applicable at strong stress gradients in the vicinity of cracks, 
or more generally, in highly heterogeneous media. We must also warn the reader: the indeterminate couple 
stress model is, in our view, a certain singular limit of the Cosserat model with independent displacements and 
micro-rotation and therefore some degenerate behaviour is to be expected throughout. 

^ Le Roux m seems to give for the first time a second gradient theory in linear elasticity using a variational formulation |65II67| . 
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1.2 The linear indeterminate couple stress model 

As hinted at above, the indeterminate couple stress model is a specific gradient elastic model in which the 
higher order interaction is restricted to the continuum rotation skew Vu (or equivalently, curlw). It is therefore 
traditionally interpreted to include interactions of rotating particles and it is possible to prescribe boundary 
conditions of rotational type. Superficially, this is the simplest possible generalization of linear elasticity in 
order to include the gradient of the local continuum rotation as a source of stress and strain energy. In this 
paper, we limit our analysis to linear isotropic materials and only to the second gradient!^ of the displacement 

D'^u = Uk,ij = Sik,j - ejk,i - where e = sym Vu. 

In general, the strain gradient models have the great advantage of simplicity and physical transparency since 
there are no new independent degree of freedoms introduced which would require interpretation. Since in this 
model there are no additional degrees of freedom (as compared to the Cosserat or micromorphic approach) the 
higher derivatives introduce a “latent-microstructure” (constrained microstructure). However, this apparent 
simplicity has to be payed with much more complicated traction boundary conditions, as will be seen later. 

We will see in Section |4l surprisingly, that the mentioned rotational interaction can equivalently be viewed 
as a strain type interaction in the indeterminate couple stress model. Therefore, the first interpretation of 
rotational interaction (which is classical) is ambiguous as long as the problem is not specified together with 
boundary conditions appearing as effect of the kind of partial integration which is performed. We may choose, 
contrarily to our intuition, another representation of the curvature energy motivated by formal considerations 
of invariance properties. In this regard we highlight the fact that force stresses for a material of higher order 
are far from being uniquely defined: it is always possible to add a self-equilibrated (divergence-free tensor field) 
force field changing the constitutive stress tensor but leaving unaltered the equilibrium equations [26l [73]. 

Often, such kind of models introduce too many additional parameters (or too many additional artificial 
degrees of freedom) which are neither easily interpreted, nor easily to be determined from experiments. Our 
discussion may also be interpreted with the background to only include those higher order terms that are 
required to describe the pertinent physics. It is clear that higher order models should not be more complicated 
than is warranted by experimental observation. A permanent nuisance in this respect is the question of how 
to identify new material parameters which are connected to the possible non-symmetry of the total force-stress 
tensor having the same dimensions as the classical shear modulus /r [N/mm^j. In the Cosserat model the 
connecting parameter is the Cosserat-couple modulus EH Il2], which, for the indeterminate couple stress 
model considered here, is formally p,c —t c». 

The Cauchy-Boltzmann axiom, well known from classical elasticity, requires the symmetry of the force stress 
tensor and may serve us also in the realm of this higher order theory to restrict the bewildering possibilities. 
Already Cauchy wrote [HI p. 344-345]: 

"... les composantes A,T,£;T,B,'D;£,'D,C des pressions supportees au point P par trois plans 
paralleles aux plans coordonnes des yz, des zx et des xy, pourront etre generalement considerees 
comme des fonctions lineaires des deplacements 77, C, et des leurs derivees des divers ordres. ’d. 

Truesdell and Toupin |1051 p. 390] write: 

““Theories of elastic materials of grade 2 or higher had been proposed by several authors [Cauchy 
m, St. Venant Jaramillo mh but under the assumption that the [total-force] stress tensor is 
symmetric”. 

Indeed, Jaramillo |47] considers a second gradient elastic material and obtains the dynamic equations by Hamil¬ 
ton’s principle. He observes dispersion relations in wave propagation problems. For simplicity only he restricts 
his discussion to those second gradient formulations, which give rise to a symmetric total force-stress tensor 
and obtains a classification for isotropic materials EZl P- 51, Eq. (96)]. The subject was pushed forward in 
the late 1950’s with works of Toupin |102[ 1103) . Grioli [SIl [35] , Mindlin and Koiter [^, among others, 
see the references later in this paper. Yang et al. [106] give an erroneous motivation for a symmetric moment 

^ There is such a formula, which says that all second derivatives of u can be obtained from linear combinations of partial 
derivatives of strain, i.e. D^u = Lin(Vsym Vu), Uk,ij = £ik,j — ^jk,i — ^ij,k, where e = sym Vu. 

^Our translation: The components [of the symmetric total force stress tensor] A, J^, £; £,'D,C can be considered in 

general as linear functions of fho displacement and their derivatives of arbitrary order. 
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stress tensor, as will be shown in [^. Neff et al. [85] considered the singular stiffening behaviour for arbitrary 
small samples in the Cosserat and indeterminate couple stress model and concluded that in order to avoid these 
singular effects one has to take a symmetric moment stress, thus providing the first rational argument in favour 
of symmetric moment stresses. In [85] the same modefl has been derived based on a homogenization procedure 
and a novel invariance requirement introduced by Neff et al. |48) . called micro-randomness and it has been 
shown that the model is well-posed. 

1.3 Our perspective 

Our contribution is intended to clarify and delineate under what boundary conditions we may expect or use 
symmetric nonlocal force stresses in the indeterminate couple stress model. When trying to relax the 4.th order 
problem (from gradient elasticity), it also seems expedient to retain the symmetry of the force stress tensor 
and of the moment stress tensor. Respecting symmetry restricts the possibilities to choose among 2.nd order 
micromorphic models. The importance of switching to a 2.nd order problem with new independent degrees of 
freedom is clear from the implementational point of view with finite elements: a 2.nd order problem is much 
easier and more efficient. However, given the antisymmetric classical and our new symmetric formulations we 
may arrive at completely different 2.nd order formulations in case of mixed displacement-traction boundary 
conditions. 

In general, the hyperstress-tensor (couple stresses, sometimes called double-stress [65]) in second gradient 
elasticity 11021 l69l EU 11001 fTHl flT] (see also the recent papers [SO] ET] [191 11071 ES] E] [TIJ EU EU EH EH ES] ) may 
be defined as m = (triyfc) = Djj2^W{D'^u). Since = (uijk) is a third order tensor, so is m. Moreover, since 
Uijk is symmetric in (jk) the same is usually assumed for triyfc- This, however, is not mandatory, see [75] . 

In the framework considered in this paper, the hyperstress-tensor is defined as 

7Ti .— ZlyITcurv (^(curl u)) Or Tn .— T^Curi (sym Vii) i^curv(Corl (sym Vzi)), 

respectively, and both expressions are 2.nd order tensor^ and are also called couple stress tensors, since they act 
as dual objects to gradients of rotations. On the other hand, as we will see, we have two competing expressions 
of the nonlocal force stress tensor: a symmetric tensor r versus an anti-symmetric tensor r: 

t = /jLI symCurl {2 oi dev symCurl (sym Vu) -I-2 02 skew Curl (symVu)} = sym Curl (m) G Sym(3) 

T = ^ antiDiv{2ai devsym V[axl(skew Vu)] -I- 2a2 skew V[axl(skew Vu)]} = ^ antiDiv(m) G so(3), 

with Div(r — r) = 0. Since [Curl(sym Vu)]^ = V[axl(skew Vu)], it follows 

r —r = /iL^| 2 (sym Curl — - antiDiv).[devsymCurl (symVu)] 

-|-2q! 2 (symCurl -fi antiDiv).[skew Curl (symVu)] |. 

The independent constitutive variable k := V(curlu) is the second gradient contribution considered by Grioli 
[37] . Toupin [102] . Mindlin [72], Koiter [51] and Sokolowski |100j . In general, neither fh nor m couple stress 
tensors are symmetric. 

The symmetry of the force stress tensor in continuum mechanics is regulary discussed in the literature, see 
e.g. ElEiES]. It has been suggested by McLennan [55] that a symmetric force stress tensor can always be 
constructed by adding divergence-free couple stresses, since only its divergence occurs in the local conservation 
law. However, all of the previously given expositions use anti-symmetric nonlocal force-stresses. Since there 

It must be noted that the grandmaster Koiter |51l p. 17-19, 23, 41] came to reject the significant presence of couple stresses 
because he based his investigations on the indeterminate couple stress theory with uniformly pointwise positive definite curvature 
energy, which tends to maximize the influence of length scale effects in its rotational formulation. His arguments only show that this 
special constrained gradient theory together with it’s boundary conditions cannot be based on experimental evidence. However, 
the main thrust of his comments remains valid and our symmetric formulation may compare favorable. We should also have in 
mind that Mindlin ceased to use these models because he could finally not see the physical relevance at this time. Truesdell and 
Noll also wrote [Tol P- 400]: “In favour of the Grioli-Toupin theory, in which the microrotation and macrorotation coincide, we 
can find no experimental evidence or theoretical advantage.” 

^See the Appendix for the relation between the second order tensor fh and the third order tensor m. 
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is no conclusive evidence for the real need of a non-symmetric total force-stress tensor in the purely mechan¬ 
ical context, we apply Ockham’s razor and discard these non-symmetric force stress formulations. Our new 
alternative formulation will have symmetric couple stresses and symmetric force stresses. Thus it satisfies the 
Cauchy-Boltzmann’s axiom. We also show that the new formulation is well-posed in statics. While conceptually 
very pleasing, the real merits of such a “completely symmetric” formulation have yet to be discovered. 

Similarly to the classical indeterminate couple stress model which can be obtained as a constrained Cosserat 
model, our new Curl (sym Vu)-model can be obtained as a constrained “microstrain” model [Ml IMl [78] . 

The question of boundary conditions in higher gradient elasticity models has been a subject of constant 
attention. Bleustein has formulated the conclusive answer for general gradient elastic models involving the 
surface divergence operator jlOj . However, the traction boundary conditions obtained by Tiersten and Bleustein 
in m with respect to the special case of the indeterminate couple stress model are incomplete. In a forthcoming 
paper [61] we discuss and correct the form of the traction boundary conditions considered until now in the 
classical indeterminate couple stress model [72l 11021 [STJ [HU |93l lU 11061 IM] . Here, we just provide the correct 
answer obtained there in the form of a summarizing box. 

The plan of the paper is now as follows: after a subsection fixing the notation, we outline some related models 
in isotropic second gradient elasticity; we prove some auxiliary results and we discuss the invariance properties 
of the considered energy; we recall the classical indeterminate couple stress model with skew-symmetric nonlocal 
force-stress (i.e. with non symmetric total force-stress tensor); we formulate the equilibrium problem for the 
new isotropic gradient elasticity model with symmetric nonlocal force stress (i.e. with symmetric total force- 
stress tensor) and we give an existence result; we discuss the difference of the classical indeterminate couple 
stress model with the introduced symmetric model; paying particular attention to the boundary virtual work 
principle we show that these two possible formulations are applicable for different types of traction boundary 
conditions; we discuss the possibility to switch from a 4.th-order problem to a 2.nd order micromorphic model. 
All our existence results can be extended, mutatis mutandis, to first order anisotropic behaviour [saiMiiiQ], 
i.e. considering as total energy (C. sym Vu, sym Vm) -I- Wcurv{D‘^u) as long as C is a uniformly positive definite 
tensor. We hnish with some boxes summarizing our models and findings. 

1.4 Notational agreements 

In this paper, we denote by the set of real 3x3 second order tensors, written with capital letters. 

For a,b G we let (a, &)r 3 denote the scalar product on R^ with associated vector norm ||a||g 3 = (a,a)R 3 . 
The standard Euclidean scalar product on R^^^ is given by (A, F)r 3 x 3 = tr(AF^), and thus the Frobenius 
tensor norm is ||A|p = (A,A)r 3 x 3 . In the following we omit the index R^,R^^^. The identity tensor on 
r3x3 ijg denoted by 1, so that tr(A) = (A, 1). We adopt the usual abbreviations of Lie-algebra theory, 
i.e., so(3) := {A G R^^^ |A^ = —A} is the Lie-algebra of skew symmetric tensors and s[(3) := {A G 
r3x3 I tr(A) = 0} is the Lie-algebra of traceless tensors. For all A G R^^^ we set sym A = ^(A^ -|- A) G Sym, 
skew A = i(A —A^) G so (3) and the deviatoric part dev A = A — ^ tr(A)l G s[(3) and we have the orthogonal 
Cartan-decomposition of the Lie-algebra g[(3) 

g[(3) = {sl(3) n Sym(3)} © so(3) 0 R-1, A = devsym A + skew A + i tr(A)l. (1.5) 

O 

Throughout this paper (when we do not specify else) Latin subscripts take the values 1, 2,3. Typical conventions 
for differential operations are implied such as comma followed by a subscript to denote the partial derivative 
with respect to the corresponding cartesian coordinate. We also use the Einstein notation of the sum over 
repeated indices if not differently specified. Here, for 

( 0 -03 02 \ 

03 0 —oi j G so(3) (1.6) 

-02 Oi 0 / 

we consider the operators axl : so(3) —>■ R^ and anti : R^ —so(3) through 

axl(A) := (oi, 02 , 03 )^ , A. u = (axl A) x u, (anti(n))y = —VuGR^, (1.7) 

(axl A )/2 = — — 2 : -^ij — (^xl A)/j; =: anti(axl A)ij, 
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where is the totally antisymmetric third order permutation tensor. We recall that for a third order tensor 
E and X G f G we have the contraction operations E : X G E. n G and X.v G with the 

components 

(E : X)i = Eyfc Xfej , (E. = Eyfe nfe , {X.v)i = XijVj. (1.8) 

For multiplication of two matrices we will not use other specific notations. 

We consider a body which occupies a bounded open set of the three-dimensional Euclidian space and 
assume that its boundary dfl is a piecewise smooth surface. An elastic material fills the domain 17 C R^ and we 
refer the motion of the body to rectangular axes Oxi. By C“(17) we denote the set of infinitely differentiable 
functions with compact support in 17. In order to realize certain boundary conditions on an open subset F C 917 
we make use of the space of functions that vanish in a neighborhood of F, i.e. 

C^{n,T):={u\3vGC^{m.’^\T) such that r'lp = u} ■ 

Here, is a vector tangential to the surface 917 \ F and which is orthogonal to its boundary 9(917 \ F) = 9r, 



Figure 1: The domain ft C R® together with the part F C dfl, 
where (geometric) Dirichlet boundary conditions are prescribed. 
We need to represent the boundary conditions on a disjoint union 
of dQ = (dfl \ r) U r U dr, where F is a open subset of 


r~ = n X h'~ is the tangent to the curve dT with respect to the orientation on d^l \ T. Similarly, i' := 
is a vector tangential to the surface F and which is orthogonal to its boundary 9r, t := = n x is the 

tangent to the curve 9r with respect to the orientation on F. The jump across the joining curve 9F is defined 
by [■ ]^ — [• ]“, where 


lim _ 

X g 9Q \ r 
X ^ sr 


[.]+:= lim [.]. 

X g r 
X -j. ar 


We assume that 917 is a smooth surface. Hence, there are no singularities of the boundary and the jump 
[•]■'■ — [•]“ arises only as consequence of possible discontinuities of the corresponding quantities which follows 
from the prescribed boundary conditions on F and 917 \ F. 

The usual Lebesgue spaces of square integrable functions, vector or tensor fields on 17 with values in R, R^ 
or R^^^, respectively will be denoted by L^(17). Moreover, we introduce the standard Sobolev spaces [T1 IB51H7] 

Hi(17) = {u G L2(17) I gradw G T^(17)}, 11^11^1(0) := l|w|li 2 (o) + IlgradulH^^j.^) , 

H(curl; 17) = {u G \ curlw G L^{n)}, |lw||H(curi;n) — Il^lli2(n) + l|curlw|li 2 (o) , (1.9) 

H(div;17) = {u G L2(17) |divu G T^(17)}, IMlH(div;n) '= Mhia) + , 

of functions u or vector fields v, respectively. Furthermore, we introduce their closed subspaces FfQ(17), 
Ho(curl; 17) as completion under the respective graph norms of the scalar valued space (^“(ll). We also consider 
the spaces 

iji(div;17;r), i?i(curl; 17; F) 

as completion under the respective graph norms of the scalar-valued space of the scalar-values space (7“ (17, F). 
Therefore, these spaces generalize the homogeneous Dirichlet boundary conditions: 

m|p = 0, and (it, n)|r = 0 and u x n|r = 0, 
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respectively. For vector fields v with components in H^(n), i.e. v = (r’l,'C 2 ,,Vi £ H^(f2), we define V v = 
((V vi)'^, (V f 2 )^, (V 1 ^ 3 )^) , while for tensor fields P with rows in H(curl; fl), respectively H(div ; fl), i.e. P = 
[Pi , P 2 , P 3 ), Pi S H(curl;r2) respectively Pi S H(div ; n) we define Curl P = ((curlPi)^, (curlP 2 )^, (curlPa)^)^ 
DivP = (divPi,divP 2 ,divP 3 )^. The corresponding Sobolev-spaces will be denoted by 

H\n), pi(Div;n), H\Cml;n), pi(Div; F), pi(Curl; F). 


2 Preliminaries 


2.1 Related models in isotropic second gradient elasticity 

One aim of this paper is to propose a new representation of the curvature energy Wcurv(P^n) and to prove that 
the corresponding minimization problem 


I{u) = [ 
Jn 


^ IlsymVwlP + - [tr(symVu)]^ + IFcurv(P^u) 


dV 




min. w.r.t. u, 


( 2 . 1 ) 


admit unique minimizers under some appropriate boundary condition. Here A, fi are the usual Lame constitutive 
coefficients of isotropic linear elasticity, which is fundamental to small deformation gradient elasticity. If the 
curvature energy has the form Wcurv{D^u) = Wcurv(Psym Vu), the model is called a strain gradient model. 
We define the third order hyperstress as Do'^uWcurv{D^u). 

In the following we outline some curvature energies already proposed in different isotropic second gradient 
elasticity models: 

• Mindlin [511 EHl El] considered energies (gradient elastic) based on the tensors 


'Hijk — 

'Hijk — Uj,ki) — ^kj,ii 

kij — ~^jlk‘Uk,li 
S 1/ 

Vijk ~ ^\'^k,ij + Uijk + Uj^ki)- 


(i.e. 77 = V(Vu)), 

(i.e. ?7 = V(symVu))), 

(i.e. k = iv(curlu)), 


The most general isotropic curvature energy defined in terms of D^u has 5 material constants, while the 
anisotropic representation is much more involved and still subject of ongoing research [7||23]. 

Mindlin and Eshel m have also proposed the following three alternative forms : 

IFcurv(P u) = /r ^ 'l^kiiVkjj T ^ "kfiji^Tlij]^ + Ug ^ klijf^T]ji^i + ^ k^jjiTji^i^i + ^ (I) 

7-2 r ( 2 )~ ~ , ( 2 )~ ~ , ( 2 )~ ~ , (2)~ ~ , (2)~ ~ 1 

= [a( 'rjuk'nkjj + a '2 ’nijilikk + % Vuklljjk + 04 ’lyfc’Iyfc + as Vijkflkji] (II) 

= fiLl + d^s'^vfijVkkj + af'^vfjkvfjk + (^ijkhjVkii], (HI) 

which are frequently cited in the literature, where Pc is the smallest characteristic length in the body and 
(7) '^(i) 

a] ,al are dimensionless weighting parameters. 

• a simple curvature energy is considered by Lam (53] [88] 


IFcurv(P>^u) = MPc[ao ||VdivM||^ + ai r]ijkr]ijk + 02 || sym V(curlu)||^] 

= /iPc[ao ||V tr(symVu)|p + m %jk%jk + 4 a 2 || symCurl (symVM)|p] 

where %jk is called the deviatoric stretch gradient and which is defined (see e.g. [85]) by 

„( 0 ) — IfA.. -I- A., r?'® A- d, - 1 

> 'lijk — C- 'immk ^ ^jk 'Immi w ivfcj 'ImmjJ- 


( 2 . 2 ) 


Vijk = vfjk - V^fk 




• another simplified strain gradient elasticity model is proposed in [Sj EU [56] based on the curvature energy 


= fiLl [ao ||Vtr(symV'u)|p + oi ||V (symVM)|p] (2.3) 

= fiLl [ao ||V divujl^ + ai jjV (symVM)||% 


which already leads to symmetric nonlocal force-stresses, see Section |3l 

• in the same line, using also the second order curvature tensor k = iVcurlrt, in [1081 [5(1] the following 
energy is considered 

Wcurv{D^u) = fxLl [ao ||Vdivu||^ -I- oi [[Vcurl u|[^] = [oo ||Vtr(sym Vu)||^ -I- oi || V[axl(skew Vu)]||^] 
= [ao ||Vtr(sym Vu)||^ -I-ai || V (skewV'u)|p]. (2.4) 


Let us remark that tr(A:) = tr(V axlskewVu) = div(curl u) = 0. 


the indeterminate couple stress model (Grioli-Koiter-Mindlin-Toupin model) [ST] [21 ETJ [72] 11031 
11001 139] in which the higher derivatives (apparently) appear only through derivatives of the infinitesimal 
continuum rotation curlu. Hence, the curvature energy has the equivalent forms 


Wcurv(-D^M) = mLc ^11 sym V(curl u) 


Ql2 

-I- — II skew V(curl u) 


= [oi II sym V[axl(skew Vu)]|p -I- 02 || skew V[axl(skew Vm)]||" 




(y.\ ‘ 

— II devsym V(curl w)||^ + ^ || skew V(curl u 


0-2 


L 4 " 4 

= [oi II sym Curl (sym Vu)|p -I- 02 || skew Curl (sym Vu)|p] 

= [oi II dev sym Curl (sym Vu)|p -I- 02 || skew Curl (sym Vu)|p] 


(2.5) 


Note carefully that tr[sym Curl (sym Vm)] = tr[sym V[axl(skew Vm)]] = 0. Therefore, we are entitled 
to use the deviatoric-representation, which is useful when regarding the model in the larger context of 
micromorphic models. Here, we have used the master identity to be established in Corollary 12.21 

V[axl(skew Vu)] = [Curl (sym Vm)]^, 

'-V-' '-V-' 

rotation gradient strain gradient 


which allows us easily to switch from considerations on the level of strain gradients to the level of rotational 
gradients and vice versa. 

We also used the identities 


2 axl(skewVu) = curia, symV(curlM) = 2 sym Curl (sym Vm), 
skewV(curlu) = —2 skew Curl (sym Vm), tr[Curl (sym Vu)] = 0. 

Although this energy admits the equivalent forms (12.51) i and (12.51) fi. the equations and the boundary value 
problem of the indeterminate couple stress model is usually formulated only using the form (12.51) i of 
the energy. Hence, we may reformulate the main aim of the present paper: to formulate the boundary 
value problem for the indeterminate couple stress model using the alternative form (I2.5l) fi of 
the energy of the Crioli-Koiter-Mindlin-Toupin model. We also remark that the spherical part of the 
couple stress tensor remains indeterminate since tr(V(curlu)) = div(curlM) = 0. In order to prove 
the pointwise uniform positive definiteness it is assumed, following m, that oi > 0 ,02 > 0. Note that 
pointwise uniform positivity is often assumed when deriving analytical solutions for simple boundary value 
problems because it allows to invert the couple stress-curvature relation. We will see subsequently, that 
pointwise positive definiteness is not necessary for well-posedness. 
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• In this setting, Grioli [37l[39] (see also Fleck [30l [HI |32] ) initially considered only the choice oi = a-i. In 
fact, the energy originally proposed by Grioli [37] is 






^ ||V(curl m)|P + ^ tr[(V(curl u))^] 

[ai II devsym V[axl(skew Vu)]|p + oi || skew V[axl(skew Vm)]|P 
+ 77 ' (V[axl(skew Vm)], (V[axl(skew Vu)])'^)] 


( 2 . 6 ) 






^ II devsym V(curl u)|p + ^ || skew V(curl u)|p + ^ (V(curl u), (V(curl u))’^) 

CH + , mi2 , 0^1 “ V 


dev sym V (curl u) 


II skew V (curl u) 


Mindlin [721 p. 425] explained the relations between Toupin’s constitutive equations [102] and Grioli’s [37] 
constitutive equations and concluded that the obtained equations in the linearized theory are identical, 
since the extra constitutive parameter rj' of Grioli’s model does not explicitly appear in the equations 
of motion but enters only the boundary conditions. The same extra constitutive coefficient appears in 
Mindlin and Eshel’s (III) and Grioli’s version (|2.6I1 . 

• the modified - symmetric couple stress model - the conformal model. On the other hand, in the 
conformal case [S3 |S1] one may consider that 02 = 0 , which makes the second order couple stress tensor 
m symmetric and trace free m- This conformal curvature case has been considered by Neff in [SS]: the 
curvature energy having the form 

WcurviD'^u) = ^ II sym V(curl m)||^ = fiL^ai || devsymGurl (sym Vm)|P. (2.7) 

Indeed, there are two major reasons uncovered in [85] for using the modified couple stress model. First, in 
order to avoid singular stiffening behaviour for smaller and smaller samples in bending [83] one has to take 
a 2 = 0. Second, based on a homogenization procedure invoking an intuitively appealing natural “micro¬ 
randomness” assumption (a strong statement of microstructural isotropy) requires conformal invariance, 
which is again equivalent to 02 = 0. Such a model is still well-posed [35], leading to existence and 
uniqueness results with only one additional material length scale parameter, while it is not pointwise 
uniformly positive definite. 

• the skew-symmetric couple stress model. Hadjesfandiari and Dargush strongly advocate [42ll43ll44] 
the opposite extreme case, Oi = 0 and 02 > 0, i.e. they propose the curvature energy 

WcurviD^u) = \\ skew V(curlu)||^ = axlskew V(curlM)||^ (2.8) 

= ||curl (curlu)lp = 02 II skewGurl (symV'u)lp. 
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In that model the nonlocal force stresses and the couple stresses are both assumed to be skew-symmetric. 
Their reasoning, based in fact on an incomplete understanding of boundary conditions (see [61j l is critically 
discussed and generally refuted in [87], while mathematically it is also well-posed. 


2.2 Auxiliary results 

Further on, we consider a simply connected domain 17 C The starting point is given by the well-known 

Nye’s formula [SU] US] 

-Gurl !■= (Vaxl3')'^-tr[(Vaxl A)"^]!, 

K := V(axl A) = —(Gurl A)^ -|- i tr[(Curl A)^] 1 = ^ tr[a] 1 “Nye’s curvature tensor”, (2.9) 

for all skew-symmetric matrices A G so(3), where a := — Curl A is the micro-dislocation density tensor. 
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Proposition 2.1. Let p : C —>• be given. The formula 

V[axlskewp] = [Curl (sym p)]^ (2-10) 

holds true if and only if there is u € C^(Ll) such that p = Vit. 

Proof. Let us first prove that 

V[axl(skew Vu)] = [Curl (sym Vm)]^, for all uGC^{LI). (2-11) 

On the one hand, using Nye’s formula for A = skew Vu, we obtain 

— Curl (skewVu) = (V[axl(skewV m)])^ — tr[(V[axl(skew Vu)])"^] 1, (2-12) 


which implies 

— Curl (skew Vu) 


(V[axl(skew Vu)])”^ — ^ tr[(V(curlM))’^] 1 

(V[axl(skew Vu)])^ — idiv(curlu) 1 = (V[axl(skew Vm)])^. 


(2.13) 


On the other hand, Curl (Vu) = 0, Vu = sym Vu + skew Vu. Thus, we deduce 

Curl (sym Vu) = Curl (Vu — skew Vu) = Curl (Vu) — Curl (skew Vu) (2-14) 

Curl(Vu) + (V[axl(skew Vu)])^ — idiv(curlu) 1 = (V[axl(skew Vu)])^. 

This establishes the first part of the claim. 

Now, we prove that V[axlskewp] = [Curl (sym p)]^ implies that there is a function u £ (7^(0) such that 
p = Vu. Using again Nye’s formula, we obtain 

(V axlskewp)^ —(Curl (skewp)) + tr[V[axlskewp]] 1. (2.15) 

Hence, our new hypothesis is Curl (symp)=(V[axl skewp])^, which implies Curl (symp)= — (Curl skewp) + 
tr[V[axlskewp]] 1. Hence, we obtain 

Curl (symp + skewp) = tr[V[axlskewp]] 1 Curl (p) = tr[V[axlskewp]] 1, (2.16) 

or, in the equivalent form 


Curl (p) = div(axlskewp) 1. 


(2.17) 


We have obtained the formula 


tr[Curl p] = 3 div(axlskewp). (2-18) 

Let us also remark that considering a matrix B £ we have 

tr[Curl (skewH)] = 2(6i^i + & 2.2 + 53 . 3 ) = 2 div 6 , & = axl(skew 5). (2-19) 

Therefore, from (12.1911 we also have obtained 

tr[Curl (skewp)] = 2 div [axl(skewp)]. (2.20) 

Moreover, for a matrix B £ R^^^, we have that 

tr(Curl B) = {Bi3^2 — 1131 , 2 ) + (ll21,3 “ lll2,3) + (ll32,l — ll23,l). 

We deduce tr(Curl S') = 0 for all S £ Sym(3). Hence, 

tr[Curl (symp)] =0 Vp £ R^^^. (2-21) 

The relations (12.201) and (12.211) lead to tr[Curl p] = 2 div [axl(skewp)], and together with (I2.18P to 

div [axl(skewp)] = 0. (2.22) 

Using (I2.17|) . we obtain Curl p = 0. Since H is an open domain in R^, it follows that there is an vector u, such 
that p = Vu and the proof is complete. □ 
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Corollary 2.2. For u G (7^(0) the following formula holds true 

V[axl(skew Vm)] = [Curl(sym Vu)]^. (2.23) 

Corollary 2.3. For u G (7^(0) the following formula holds true 

[Vcurlw]^ = Curl [(Vu)^]. (2.24) 

Therefore, (Vm)^ G i7(Curl;f2) is equivalent to curlw G 77^(12). 

As consequence of the above remark, it follows that if curlu G i7^(n), then {Vu)’^.t G Lf{dVl) for any 
tengential direction r at the boundary and, since {{Vu)'^ .T,n) = (t, {Vu).n), it results that (l — ntSm) {Vu).n G 
Lf{dV>.), in the sense of trace. 

Let us also recall the Saint-Venant compatibility condition 
Proposition 2.4. (see e.g. [14]) Let a symmetric tensor field e : 12 C —)> Sym(3) be given. Then, 

inc(e) := Curl [(Curl e)^] =0 O there is uG (7^(12) such that e = sym(VM). (2.25) 

We note that 

inc(p) := Curl [(Curl symp)^] G Sym(3) (2.26) 

while Curl [(Curl skewp)^] e so(3) V p S 

We also remark that a direct consequence of Proposition 12 .1 1 is the following first order compatibility condition 
Proposition 2.5. Let p : 12 C R^ —>• R^^^ be given. Then, 

INC(p) := [Curl (sym p)]^ — V[axlskewp] there is uGC‘^{LI) such that p = Vu. (2.27) 

We observe that 

INC(p) e R3^3 y ]^3x3_ 

We recall the well known first order compatibility condition 
Proposition 2.6. Let p : 12 C R^ —>■ R3x 3 jg giygn. Then, 

Curl p = 0 ([Curl (symp)]^ =—[Curl (skewp)]^) there is m S (7^(12) such that p = Vu. (2.28) 

Hence, we have the following equivalence 
Corollary 2.7. Let p : 12 C R3 —R3 x 3 ^g Then, 

Curl p = 0 INC(p) = 0. (2.29) 

We finally remark that 

Curl(INC(p)) = Curl [(Curl symp)^] = inc(symp). (2.30) 

2.3 Discussion of invariance properties 

The difference between the V[axl(skewVu)] formulation and the Curl(symVu) formulation can be seen when 
considering the results under superposed incompatible tensor fields: 

Remark 2.8. The quantity [Curl (sym Vit)]^ is invariant under locally adding a skew-symmetric non-constant 
tensor field W(x) G so(3), i.e. 

[Curl (sym(Vu-I-W(a:)))]^ = [Curl (sym Vw)]^. (2-31) 

However, since V[axlskewp] ^ [Curl (symp)]^ for general incompatible p G p Vu, the quantity 

V[axl(skew Vw)] is not invariant under locally adding W{x) G so(3), i.e. 

V[axl(skew(Vu -I- W{x))\ ^ V[axl(skew Vu)]. (2.32) 
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Remark 2.9. The term V[axl(skew Vu)] is invariant under locally adding a symmetric, non-constant tensor 
field S{x) £ Syin(3), i.e. 


V[axl (skew(Vu + 5'(a:))] = V[axl(skew Vm)]. (2.33) 

Let us recall the Lie-group decomposition GL~*’(3) and the corresponding Lie-algebra decomposition: 

GL'''(3) = {SL(3)/SO(3)} • SO(3) • (R'*"-!) Lie-group decomposition, 

TiGL'''(3) = = fll(3) = {s[(3) fl Sym(3)} © so(3) © R-1 Lie-algebra decomposition. (2.34) 

The space Sym(3) is not a Lie-algebra, it is only a vector space and it does not have a group structure: the set 
GL(3)/SO(3) = PSym(3) is not a group, neither is the set fll(3) fl Sym(3) a Lie-algebra. Hence, the invariance 
requirement in (12.311) . i.e., locally adding W{x) G so(3) is much more plausible than assuming (12.331) since it 
yields so(3)-Lie invariance. 

2.4 Conformal invariance of the cnrvature energy and group theoretic arguments 
in favour of the modified couple stress theory 

An infinitesimal conformal mapping [821185] preserves (to first order) angles and shapes of infinitesimal figures. 
The included inhomogeneity is therefore only a global feature of the mapping (see Figure [2]). There is locally 
no shear-type deformation. Therefore it seems natural to require that the second gradient model should not 
ascribe energy to such deformation modes. 

A map (/)c : R^ —!> R^ is infinitesimal conformal if and only if its Jacobian satisfies pointwise Vfcix) G 
R- 1 +so(3), where R- 1 +so(3) is the conformal Lie-algebra. This implies [511|5S1[H3] the representation (see 
Figure [2]) 

4>c{x) = - (2(axl VF, x) x — axlkF||x|p) + [pi A\ - x -\-b, (2.35) 

where W,A G so(3), b G R^, p G R are arbitrary given constants. For the infinitesimal conformal mapping fc 
we note 

V(j)c{x) = [(axl W, x) + p] 1 + anti(JF. x) + A, 
skewV(/)c(x) = anti(kF.x) + A, 
devsym V(()c(x) = 0, 
sym Vcurl (()c(x) = 0, 

These relations are easily established. By conformal 
the curvature energy vanishes on infinitesimal conformal mappings. This is equivalent to 

Wcnrv[D‘^4’c) = 0 for all conformal maps fc (2.37) 

or in terms of the second order couple stress tensor m := i4vcuriMfFcurv(Vcurlu), 

= 0 for all conformal maps fc- (2.38) 

The classical linear elastic energy still ascribes energy to such a deformation mode, but only related to the bulk 
modulus, i.e., 


div0c(x) = tr[V0c(x)] = 3[(axlVF,x) +p], 
symV0c(x) = [(axlfF,x)1, 

Vcurl(()c(x) = 2 JF G so(3), 
skew Vcurl(/)c(x) = 2W. 

invariance of the curvature energy term we mean that 


Riin(V^c) = pII devsymy^elP© ^^ ^ [tr(V0c)]^ = ^ [tr(V0c)]^- (2.39) 

^ \ ^ z z 

^0 

In case of a classical infinitesimal perfect plasticity formulation with von Mises deviatoric flow rule, conformal 
mappings are precisely those inhomogeneous mappings, that do not lead to plastic flow m, since the deviatoric 
stresses remain zero: devsymV(()c = 0. 
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Figure 2: Infinitesimal conformal mappings [83| preserve locally angles and shapes. The corresponding couple stress tensor m is 
zero for these deformation modes. 


In that perspective 

conformal mappings are ideally elastic transformations and should not lead to moment stresses. 

Using the formulas (12.361) . it can be easily remarked that ||V[devsym Vm]|P, || devsym Vm|P, 
II devsym V(curlu)|p, || symCurl (sym Vm)|P = |; || sym V(curlM)|p are conformally invariant. Let us note 
that, using Lemma l2.11 we have 

II symCurl(sym Vu)||^ = || sym V[axl(skew Vu)] |p = ^ II sym V(curlM)||^, (2.40) 

II skew[Curl (sym Vm)]||^ = || skew V[axl(skew Vm)]|P = ^ || skew V(curlu)|p. (2-41) 

Hence || symCurl (sym Vu)|p = j || sym V(curlu)|p is also conformally invariant (use (l2.36ll fiL while 

II skew[Curl(symV(()c)]|P = I || skew V(curl(()c)r = l|W^lP> (2.42) 

and therefore || Curl(sym Vu)|p is not conformally invariant, nor is ||V(axl(skew Vu))|p conformally invariant. 
Nor is IIV tr(sym Vu)|p = ||Vdivu|p conformally invariant. 

The underlying additional invariance property of the modified couple stress theory is precisely conformal 
invariance. In the modified couple stress model, these deformations are free of size-effects, while e.g. the 
Hadjesfandiari and Dargush choice would describe size-effects. In other words, the generated couple stress 
tensor m in the modified couple stress model is zero for this inhomogeneous deformation mode, while in the 
Hadjesfandiari and Dargush choice m is constant and skew-symmetriqj. 

2.5 The classical indeterminate couple stress model based on 

II V[axl(skew Vti)] IP with skew-symmetric nonlocal force-stress 

We are now re-deriving the classical equations based on the V[axl(skew VM)]-formulation of the indeterminate 
couple stress model. This part does not contain new results, see, e.g., m for further details, but is included 
for setting the stage of our new modelling approach. 

Taking free variations 6u G C°°(D) in the energy W{e,k) = Wiin(e) -I- Wcurv(fc), but using the following 
equivalent curvature energy based on A: = V[axl(skew Vu)] = iV(curlu): 

Wcurv(A:) [o^i II devsym V[axl(skew Vu)]|p -I- 02 || skew V[axl(skew Vm)]|P] (2.43) 

=fJ,Ll [ai II dev sym Curl (sym Vu)]|p -I- 02 || skew Curl (sym Vu)]|p], 

®This observation is a further development in understanding why the Hadjesfandiari and Dargush [miiiiiii] choice is rather 
meaningless, while mathematically not forbidden [87]. 
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we obtain the virtual work principle taking free variations 6u G in the energy (12.4311 


f WCVu + t'V6u)dv 
dt Jn 


t=o 


2^ (sym Vm, sym V(5u) + Atr(Vu) tr(Vi5M) 


/rL^[2ai (devsymV[axl(skewVM)],devsymV[axl(skewV5u)]) (2.44) 


dv = 0. 


+ 202 (skew V[axl(skewVM)],skew V[axl(skewV<5M)])] + {f,Su) 

The classical divergence theorem leads to 

/ (D'rv {(T— t) + f,Su) dv — / {{a — t) . 6u) dv — / (m. n, axl(skew Vi5it))da = 0, (2-45) 

Jn Jan Jan 

for all virtual displacements 6u G where n is the unit outward normal vector at the surface dfi, a is 

the symmetric local force-stress tensor 


(T = 2 /r sym Vm + A tr(Vu) 1 G Sym(3) 

and T represents the nonlocal force-stress tensor (which here is automatically skew-symmetric) 

T = [ai antiDiv(devsym V[axl(skew Vm)]) + 02 antiDiv(skewV[axl(skew Vu)]) 
= antiDiv(devsym V(curlit)) + ^ antiDiv(skew V(curlM))] 


(2.46) 


(2.47) 


= /r Lp anti Div 


0^2 ^ 

— devsym V(curlu) + — skewV(curlu) = -antiDiv[TO] G so(3), 


where 


fh = fiL^ [oi sym V(curl u) + 02 skew V(curl u)] 

= fiL^ [oi dev sym V (curl u) + 02 skew V (curl u)] 

= fiL^ [2ai devsym V[axl(skew Vu)] +2^2 skew V[axl(skew Vu)]], 


(2.48) 


is the hyperstress tensor (couple stress tensor) which may or may not be symmetric, depending on the material 
parameters. 

The non-symmetry of force stress is a constitutive assumption. Thus, if the test function Su G C'°°(n) also 
satisfies axl(skew Vi5m) = 0 on T (equivalently curlew = 0), then we obtain the equilibrium equation 


Div 2/r sym Vm + A tr(Vu) 1 

'-V-' 

symmetric local force stress (TGSym( 3 ) 

— /rLg [anti{Div(ai devsym V[axl(skew Vm)] + 02 skew V[axl(skew Vu)])}] > +/ = 0, 

'-'"r-' ^ 

hyperstress ^m€gl( 3 ) 

"-^" 

completely skew-symmetric nonlocal force stress tGso( 3 ) 


(2.49) 


or equivalently 


Div CTtotai + / = 0, (2.50) 

where 

CTtotai = CT - f ^ Sym(3). (2.51) 

The complete consistent boundary conditions for this formulation is presented for the first time in [SIIIHD] and 
recapitulated in Figure [Hand Figure [71 
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3 The new isotropic gradient elasticity model with symmetric non¬ 
local force stress and symmetric hyperstresses 


As independent constitutive variables for our novel gradient elastic model we choose now 

e = symVu, fc = Curl (symVu) = Curl e. (3.1) 

We use again the orthogonal Lie-algebra decomposition of 

Curl (sym Vm) = dev sym(Curl (sym Vu)) -I- skew Curl (sym Vu). (3.2) 

The term i tr(Curl (sym Vit)) 1 is missing since tr(Curl (sym Vu)) = 0 anyway (already tr(Curl S) = 0 for 
S G Sym(3)). The model is derived from the free energy W{e, k) = Wiin(e) -I- Wcurv(fc), with 

.X . . _ \ 19 111 _ 119 rc 


kLiin(e) = M II symVwlc -I- — [tr(symVM)] = /r || devsymVw|c -I- — [tr(symVu)]^ 
Wcurv{k) = ai II devsymCurl (sym Vu)|p-|-a 2 II skewCurl (sym Vm)||' 


(3.3) 


conformally invariant 


not conformally invariant 


_ 2 /J+3A jg infinitesimal bulk 


where ai, a 2 are non-negative constitutive curvature coefficients and k = 
modulus, while p is the classical shear modulus. 

The hyperstress-tensor (moment stress tensor, couple stress tensor) 

in := Il^lTcurv(fc) = dev sym Curl (sym Vm) -I- 2 a2 skew Curl (sym Vrt)] 

is symmetric in the conformal case q ;2 = 0, while the nonlocal force stress tensor is always symmetric, see eq. 

ms)- 

Due to isotropy, the curvature energy Wcurv(^) involves in principle only 2 additional constitutive constants. 
Taking free variations Su G C°°(f2) in the energy (13.3p . we obtain the virtual work principle 


1 /* 

— / W{Vu + t'i76u)dv 
dt Jn 


t=o 


>n L 


2/r (sym Vm, sym Vdw)-f Atr(Vu) tr(V(5u) 

- pL^[2 ai (dev sym Curl (sym Vm), dev sym Curl (sym Vdw)) 


(3.4) 


2 02 (skew Curl (sym Vw), skew Curl (symV^u))] -I- {f,Su) 


dv = 0, 


where / is the body force per unit volume. We have the formulas 

div(i^i Qi) = {Qi, V (fii) + ipi div Qi not summed, (3.5) 

div (i?j X Si) = (S'i,curl Ri) — (i?i,curl Si) not summed, 

for all -functions : D —>■ R and Qi,Pi,Si : D —>■ R^, where ipi are the components of the vector (p and 
Qi, Pi, Si are the rows of the matrix Q, P and S, respectively, where x denotes the vector product. If we take 
in (|3.5|) Ri = [dev sym Curl (sym Vu)]i, S'; = (symV(5u)i, we get 

3 3 

^ div ([dev sym Curl (sym Vu)]i x (symVdM)^) =^((symVi5M)i,curl [dev sym Curl (sym Vm)]^) (3.6) 

2—1 2—1 

3 

— ([dev sym Curl (sym Vrt)] j, curl (sym V(5u)i). 

i=l 

Hence, we obtain 
3 

div ([dev sym Curl (sym Vu)]i x (symVi5M)i) =( (sym V(5u), Curl [dev sym Curl (sym Vm)]) (3.7) 

1=1 

— ([dev sym Curl(symVM)],Curl (symVdit)). 
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Doing a similar calculus, but choosing Ri = [skew Curl (symVu)]i, 5^ = (symV<5M)i, we obtain 
3 

^ div ([skew Curl (symVu)]i x (symV(5M)i) =((sym V5u), Curl [skew Curl (symVu)]) (3.8) 

i=l 

— ([skew Curl (sym Vu)], Curl (symVJit)). 


The above formulas lead, for all variations 5u G to 

/ ai (devsymCurl (symVu), Curl (symVJM)) + a 2 (skewCurl (symVu), Curl (syinV^w)) 


dv 


ai (symCurl [devsymCurl (sym Vu)], V(5 m, ) + a 2 (symCurl [skewCurl (symVu)], VJw,) (3-9) 


-^div 


L 

-tdiv 


Qfi ([devsymCurl (symVM)]i x (symVJu)i) + 02 ([skewCurl (symVM)]i x (symV<5M)i) 
(symCurl [devsymCurl (symVu)], V(5 m) + 02 (symCurl [skewCurl (symVM)],V5u) 

oi ([devsymCurl (symVM)]i x (symV5u)i) + 02 ([skewCurl (symVu)]i x (symVJM)^) 


dv 


i=l 


dv. 


Therefore, using the divergence theorem and a special format of the partial integration which is suggested by 
the matrix Curl-operator, it follows thalQ 


oi (devsymCurl (sym Vm), Curl (sym V^u)) -I- 02 (skewCurl (sym Vw), Curl (sym V5u)) 


dv 


(3.10) 


oi (symCurl [devsymCurl (symVit)], V(5u) -I- 02 (symCurl [skewCurl (sym Vu)], Vi5m) 


dv 


Idn _ 


^(oi ([devsymCurl (symVM)]i x (symV(5u)i) -I-Q; 2 ([skewCurl(symVM)]i x (sym V(5u)i), n)da 


Jn . 

"/< 

Jdn L 


ai (Div{symCurl [devsymCurl(symVit)] -I- 02 (symCurl [skewCurl(sym Vu)]},(5ri) 


dv 


ai symCurl [devsymCurl (symVu)] - 1-02 symCurl [skewCurl (symVrt) 


. n,6u)da 


r 3 


^^(oi [dev sym Curl (symVM)]i -I- 02 [skew Curl (symVM)]i, (symV5u)i x n)da, 

■ i=l 

where n is the unit outward normal vector at the surface did. Hence, the relation (13.41) leads to 
J Div |2/r sym Vu-I-A tr(Vu) 1 


(3.11) 


^ Lg [2 oi sym Curl [dev sym Curl (sym Vu)] -1-2 02 sym Curl [skew Curl (sym Vu)]] f + f, Su) dv 


Idn 


2fj. (sym Vm -I- A tr(Vu) 1 


.n, 6u)da 


- f {‘2dLl 

Jan 


oi symCurl [devsym Curl (sym Vm)] -1-02 symCurl [skewCurl (sym Vu)] 




.n, 5u)da 


+ 


Ian . 


(oi [dev sym Curl (sym Vu)]i -I- 02 [skew Curl (sym Vu)]i, (sym V(5M)i x n)da = 0, 


^This is an extra constitutive assumption since it is finally the form of the partial integration that determines, on the one hand, 
which force-stress tensor is generated and, on the other hand, which boundary condition is obtained. It is only the Curl-operator 
that seems to suggest this choice-but it remains a choice! 
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for all variations Su G C°°{^1). 

We can write the above variational formulation, for all variations 6u G in the following forn|l 

3 

/ (Div {a+ 9) + f,du) dv — / {{a+ 9).n,du)da — / x n, {sym'VSu)i)da = 0, (3-12) 

Jfl J dQ 


where 


<7 = 2fi sym Vu + A tr(VM) 1 G Sym(3), (local force-stress) 

9 = fiL‘1 sym{2aiCurl [dev sym Curl (symVw)] -|-2 q;2 Curl [skew Curl (sym Vu)]} (3.13) 

= symCurl (m) G Sym(3), (non-local force stress) 

TO =/iL^ {2 ai devsym[Curl (symVu)]-I-2 a 2 skew[Curl (symVu)]} G 0l(3 ) (g Sym(3) for 02 = 0). 


We call a the local force stress tensor, 9 the non-local force stress tensor and fh = Il^Wcurv(fc) the hyperstress 
tensor (couple stress tensor). 

Thus, if the test function 5u G also satisfies (symVdM)^ x n = 0 (or equivalently (symVdu).r = 0 

for all tangential vectors r at T), then we obtain the equilibrium equation 


Div{ 2^ sym Vu-I-A tr(VM) 1 

'-V-' 

symmetric local force stress cr G Sym(3) 

+ symCurl [2 ai dev symCurl (sym Vu)-f 2 q ;2 skew Curl (symVu)]}-f/ = 0. 
'---' 

hyperstress m€0l(3) 


symmetric nonlocal force stress r G Sym(3) 


(3.14) 


The first impulse is to prescribe on T C dVL the following geometric boundary conditions 


u = vP on r, (3.15) 

(1 — n C) n) (sym Vu)i X n = (1 — n 0 n) (sym Vu°)i X n, *=1,2,3 on T, 

where —>■ is a prescribed function (i.e. 3-|-2-|-2-|-2=9 boundary conditions), witl(^ it° G i?^’^(f2), and 

the following traction boundary conditions on dVl \ T 

{a + 9).n = g, on 911 \ T, (3.16) 

(TL — n 0 n) fhi X n = — n'S) n)gi, * = 1,2,3 on 912 \ T, 

where g,gi : R^ —^ R^ are prescribed functions (i.e. 3-|-2-|-2-l-2=9 boundary conditions). 

However, we need to separate normal and tangential derivatives of the test function 6u in (13.121) which is 
standard in general strain gradient elasticity, since tangential derivatives of Su are not independent of 6u. Let 
us define the matrix 


M ■= 


/ TOi X n_\ 

1 j 

1' rhi -\ 

TO2 X n_ 

, where to := 


\ TO3 X n_/ 

' 1 

1 

V TO3- / 


With the help of this matrix M, we may write 


3 

X n, (symVdit)i) = (M,symVdu) = (symM,Vdit). 

i=l 


(3.17) 


(3.18) 


^{a X b, c) = —{a, c X b). 

9lt is always possible to construct a function uq taking on the desired boundary values and having the needed regularity by 
solving 


/q (II sym Vrt|p + ||D^u|P) da —> 
Up = 
curlu „ = 


min. u, 

“0i 

curl uo Ip 


u g 
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At this point, it must be considered that the tangential trace of the gradient of virtual displacement can be 
integrated by parts once again and that the surface divergence theorem can be applied to this tangential part of 
Vdu. Before doing so, one needs to introduce (see also [MUMHH] for details) two second order tensors T and Q 
which are the two projectors on the tangent plane and on the normal to the considered surface, respectively. As 
it is well known from differential geometry, such projectors actually allow to split a given vector or tensor field 
in one part projected on the plane tangent to the considered surface and one projected on the normal to such 
surface (see e.g. [21] )■ Let {t^, r^} be an orthonormal local basis of the tangent plane to the considered surface 
at point p and let n be the unit normal vector at the same point. We can introduce the quoted projectors as 

T = ® ® = 1 — n® n {Tij = 6ij — riirij), Q = n®n {Qij = rurij). (3.19) 

In our abbreviations, the surface divergence theorem means [401 p. 58, ex. 7] 


/ {T.v,v)da= / {v,T.v)da= / {v,v)ds, (3.20) 

JdS JdS JdS 

for any field n S and v = t x n. Regarding the boundary conditions, similar as in we obtain 

3 

{fhi X n, (sym Vi5u)i) = (symM, V5u 1) = (symM, V5u {T + Q)) 

*=i = ((symM), (Vdw) T) + ((symM), iySu) Q) = ((symM) T, V5u) + ((symM), {V5u) Q) 

= ((symM) TT, V5u) + ((symM), (Vdw) Q) = {T, T(symM)V(5w) + ((symM), (VSu) Q) 

= (T, r(symM)Vdu) + ((symM), {VSu) n ® n). 


The last term on the right hand side may be rewritten in the form 

((symM), {\76u)Q) = ^(M, {\76u)Q + Q{\76u)'^) = i(M Q,\7Su) + ^{M,Q{\7Su)'^) 

= -^{M n®n, VSu) + ^ ® niVSu)'^). 

Thus, we deduce 

3 

{fhi X n, {symVSu)i) = {T,T{sYmM)VSu) + ((symM), {V6u)Q) 

= (T, T(symM)Vdu) + ((symM), (VSu) n ® n) = {T, T(symM)VdM) + ({n (8 [(V(5u).n]}(symM), 1) 
= (T, T(symM)Vdu) + {n® {(symM)(VdM).n}, 1) = {T,T{symM)VSu) + {n, {sym M)(VSu).n) 

or 

= {T,T{symM)VSu) + ((symM).n, {VSu).n) = {T,T{symM)VSu) + ((symM).n, -^) ■ 

on 

We can therefore recognize in the last term of this formula that the normal derivative 


— Su = (VSu).n = 
an 


{Sui^hnhji 


(3.22) 


of the test function field Su (the virtual displacement) appears. As for the other term, it can be manipulated 
suitably integrating by parts and then using the surface divergence theorem (13.201) , so that we can finally write 
the last summand from ()3.12l) in the form 


» 3 » 3 » 3 

/ / ,{fhi X n, {symVSu)i)da = / y {fhj x n, {sym.VSu)i)da + / y {fhj x n, {sym VSu)i) da. (3.23) 

JdQ Jr Jdn\T 

We deduce by gathering the results in (I3.21I) - (I3.23I) that the last integral on the right hand side is given by 

/ yy (mi X n, (symVdu)i)da = / {T,T{sym.M)VSu)da + / {{symM).n,{VSu).n)da. (3.24) 
Jdn\T Jdn\T Jdn\T 
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Since 


{V[T{sym M)Su]}ik = {T{syinM)).Su}i^k = {{T{symM)) {Su)j}^k 


6R3 


= {T{symM))..^^{6u)j + {T{symM)) ^^{Su)j,k 
= (T(symM))^^. + {T{symM)) ..{\/6u)jk 

= {T{syvaM))^^^^{5u)j + {r(symM) S/5u}ik, 


(3.25) 


we obtain 


[ 'S^{mi X n,{sym.VSu)i)da = [ {T,V[T{symM).Su])da — [ Tik(T{symM)). {Su)jda (3.26) 
^a\r^ ^n\r '-.-' JdQ\T 

-r„ce divergence is to be used 

((sym M).n, {'V6u).n)da. 


Jan\r 

In order to write in a compact form the above relation, let us remark that 

(T(symM)) . = Tu{symM)ij = Tii{symM)ji = Tu{syu\M)ji = {symM)jiTii = {{symM)T}j 


and further that 


{T{symM))..i^Tik = {{syinM)T}ji^kTik = (V[(symM)T] : T)^.. 


We obtain 


r,fc(T(symM))^^. Jfe), = (V[(symM)T] : T)^{6u)j = (V[(symM)T] : T,Su). 


(3.27) 


(3.28) 


We deduce by gathering the results in p.21l) - p.23l) that the last integral on the right hand side is given by 


/ ^ {sy m \7 6 u)i) da 

/an\r 


/9n\r '■ 


(r,V[T(symM).du])da- / {W[{symM) T] : T,Su)da 


surface divergence is to be used 


idn\r 


/an\r 


((sym M).n, {VSu).n)da 


{\/[{sym M) T] : T,Su) da + / {{syuiM).n, {'VSu).n)da + / ([(symM)] .i/ ,6u)ds. 

an\r Jan\T Jdian\T) 

In the above computation V[(symM) T] is not a matrix, rather a third order tensor and V[(symM) T] : T G 
is a contraction operation, i.e. 

{V[(symM)T] : T}i = {(anti[(symM) 

Similar, we handle the corresponding integral on i9r from ()3.23|1 
3 


X n, {sy ui'VSu)i) da 


(3.29) 


= — / (V[(symM) T] : T, ^u)da + / ({symM).n, {VSu).n)da + / ([(symM)]''‘.i/+, (5M)(is. 


lar 
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Therefore, the variational formulation (I3.12|) can be rewritten as 


f (Div (cr + r) + /, 6u) dv — j 

J fl J 


( (cr + f). n — V[(symM) T] : T 

'-,-^ 

6u — surface contact forces depending 
on the curvature of the boundary 

independent first order variation 


, Su)da 


(3.30) 



((symM).n, (V(5rt).n) da — 

V ^ 

double force 

"-V-" 

Su — independent second order 
normal variation of the gradient 



(([symM]”*" — [symM] ).u,6u)ds = 0, 

^ ■ V' 

“edge line force” 


for all variations Su G (7°“(n), where we have used that for the regular surface dfl it holds = —v = v. 
Moreover, we also obtain 


Hence 


M.n = 



/ fhi X n S 



'/ 


7712 X 77 

.77 




y 777.3 X 77 y 



\ 


(fhi X n, n) 
(to 2 X n, n) 
{fhs X n, n) 


= 0 . 


(symM).n = 


(3.31) 


(3.32) 


On the other hand, we deduce 

.n = T .n + Q .n = (1 — n 0 n) M^.n + n 0 n .n 

= (1 — n 0 n) M'^.n + [n 0 n M'^].n = (1 — n 0 n) M^.n + n 0 [M. n].n 



/ TOl X 77 ' 



/ {fhi X 77 , 77 ) \ 

(1 — 77 0 77 ) M'^ .n + 77 0 

7772 X 77 

.77 

.77 = (1 — 77 0 77 ) M'^ .n + 77 0 

(7772 X 77 , 77 ) 


_ \ 7773 X 77 ) 



_ \ (7773 X 77 , 77 ) J 


= (1 — n 0 n) M'^.n + [n 0 0].n = (1 — n 0 n) .n . 


In view of (I3.32p . we see 


(symM).n = (1 — n 0 n) (sym M).n . 


Therefore, finally we get from p.l2l) 


(3.33) 


(DW (a+r) + f^Su) dv — / ((cr + r). n — V[(symM)T] : T, 5u) da 
2 Jdn^ - V -' 

independent first order variation 


(3.34) 


idn 


((1 — n 0 n)(symM).n, (1 — n 0 n)(V<5M).n) da — / (([symM]'*'— [symM] ).iy,Su) ds = 0, 


Su — independent second order 
normal variation of gradient 


/ar' 


“edge line forces” 


for all variations Su G An equivalent form, replacing simply r = sym Curl (m), is 

/ (Div (cr+sym Curl (m))+ /, dct) dw — / ((cr + sym Curl (m)). n — V[(sym M) T] : T, da) da 

Jn Jan' ---' 

(5u—independent first order variation 


(3.35) 


Ian' 


((1 — n 0 n)(symM).n, (1 — n 0 n)(Vdw).n) da — / (([symM]+— [symM] ).iy,Su) ds = 0. 


Su — independent second order 
normal variation of the gradient 


/ar' 


“edge line forces” 
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3.1 Formulation of the complete boundary value problem 

3.1.1 Equilibrium equation 

In terms of the symmetric force-stress tensor a and of the nonlocal force-stress tensor r which is also here 
symmetric, while the hyperstress fh G fll(3) is symmetric only for 02 = 0, the equilibrium equations may now 
be written in the formalo 


DivCTtotai + / = 0, (3.36) 

where the symmetric total force stres0 is given by Ptotai = o' + t G Sym(3). 

3.1.2 Geometric (essential) boundary conditions 

To the above equilibrium equation, we adjoin on T C dfl the following boundary conditions 

u{x) = u°{x) on r, (3 be) (3.37) 

[(1 — n (8) n)(VM).n] (x) = [(1 — n (8) n)(Vu°).n](a;) on T, (2 be) 

where —>■ is a prescribed function (i.e. 3-|-2=5 boundary conditions). We assume that G 

for simplicity and transparency. If G then all tangential and normal traces of at T exist. 

Therefore, we may evaluate Vu^ at T. 

3.1.3 Traction boundary conditions 

Corresponding to the geometric boundary conditions, we have to prescribe the following traction boundary 
conditions 

{(cr-|-f).n — V[(symM) (1 — n (8) n)] : (1 — n (8) n)} (x) = t(x), ] (3 be) 

_ > a; G \ r 

[(1 — n (8) n)(symM).n] (x) = [(1 — n (g) n) g] (x), ) (2 be) 

(3.38) 

{([symM]+— [symM]“). p} (x) = 7f(x), x G OF (3 be) 

where t, g : ^ are prescribed functions on dfl \ T (i.e. 3-|-2=5 boundary conditions), while tt : —>■ R.^ 

is prescribed on 9r = d{dfl \ T) and leads to 3 boundary conditions on ST. 

Remark 3.1. If T = dQ then the solution in the C'arl [sym'Vu)-formulation and in the V[axl(skew Vu)]- 
formulation are the same, since the Euler-Lagrange equations are the same and the geometric boundary con¬ 
ditions are the same. Differences appear only if T dQ due to different specifications of traction boundary 
conditions. 


3.2 Existence and uniqueness of the solution in the Curl (sym Vti)-formulation 

In the linear couple stress theory with constrained rotations, Hlavacek and Hlavacek [46l Remark 2, p. 426] 
recognized the couple stress model already in the form (11.41) but did not give an existence result. There are 
many existence and uniqueness results for the indeterminate couple stress model in its classical anti-symmetric 
formulation. Recently, optimal results have been obtained in [451 SS] • In this section we establish an existence 
theorem for the solution of the boundary value problem {V) defined by (I3.36|) . p.37l) and p.38|) . where t = 0, 
g = 0, h = 0, = 0 and (1 — n 0 n){VvP).n = 0 for simplicity only. 

Lemma 3.2. Let u G Hj!j(f2;r) be such that symCurl (sym Vu) G L^(r2). Then, Curl(symVu) G L^(r2) and 
there is a positive constant such that 


In 


I symVu|p -I- II symCurl(sym Vm)||" 


dv > 


In 


II sym Vm|P -I- II Curl (sym Vu)||" 


dv. 


'^'^Here, infinitesimal frame-indifference amounts to W{Vu) = W(Vu-i- W),'rfW G 30(3), which is obviously satisfied. 
^^Vidoli et al. call this tensor the “effective stress tensor” [23]. 


(3.39) 
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Proof. For u G Hj(ri;r), the first Korn’s inequality implies that there is a positive constant c+ such that 

f II sym Vitlpdt! = f [^i|| sym Vm|P + i|| sym dr; > i f || sym Vujpdn + ^ f \\\7u\\'^dv. (3.40) 

Jo, Jn\‘^ 2 J ^ Jn 2 Jf 2 


On the other hand the orthogonality of sym and skew implies 

I ||Vu||^(iz) > j ||skewVu||^df. 

Jq Jq 

Therefore, there is another positive constant c"*" such that 


(3.41) 


/a L 


sym Vu||^ + II sym Curl (sym Vu)||" 


dv = 


i|| symVull^ + i|| sym Vu||^ + || sym V[axl(skew Vu)]|| 


> c"* 


> c"* 


IQ 


IQ 


I symVulP + ||skewVu|p + || symV[axl(skew Vm)]||" 


dv 


dv 
(3.42) 


symVulP + II axlskewVulP + || sym V[axl(skew Vm)]||' 


dv. 


Moreover, since axlskewVw S L^(0), the second Korn’s ineaualit\F^ (without boundary conditions and applied 
to axl(skew Vu)) implies the existence of a positive constant c’*' such that 


/n 


II axlskewVulP + || sym V[axl(skew Vm)]||^ 


dv > 


iQ 


axlskewVulP + ||V[axl(skew Vu)]|p dv. (3.43) 


Thus, there are positive constants c'^,Ci such that 


IQ L 


sym Vm||^ + II sym Curl (sym Vm)|| 


^ dv>ct 


> c"* 


In L 


IQ 


II sym Vu||^ + II axlskewVull^ + || V[axl(skew Vm)]||^ 
II sym Vu||^ + II axlskewVwll^ + || Curl (sym Vm)||^ 


In 


II symVuf + II Curl (sym Vu) 11^ 


dv. 


The proof is complete. 


dv 
dv 

(3.44) 

□ 


Let us consider that we have considered null boundary conditions for simplicity. Hence, in the following we 
study the existence of the solution in the space 


On Xq we define the norm 


Xo = {uG iLp (H; r) I sym Vm G iL(Curl; O)}. 


Mxo = (l|VM||i 2 (o) + II Curl (sym VM)||i 2 (o) 


(3.45) 


(3.46) 


and the bilinear form 


{u,v) = [ : 

JQ . 


2 fi (sym Vu, sym Vd) + A tr(Vu) tr(Vu) 

+ /r [2 ai (dev sym Curl (sym Vw), dev sym Curl (sym Vu)) 


(3.47) 


+ 2a2 (skew Curl (sym Vu), skew Curl (symV?;))] 


dv. 


Since curlu is divergence free we also have the following Maxwell type inequality 019]; 

t| Vlcurlw] ||^2(Q) < CM(||curl (curln)||^2(Q) , for it G {it g Hq{Q) \ curl it G //(curl; fl)}. 
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where u,v € Xq. Let us define the linear operator Z : Aq —>■ R, describing the infiuence of external loads, 
i(u) = J^if, v)dv for all w G Xq. We say that w is a weak solution of the problem (V) if and only if 

{u,v) = l{v) for all v G Xq. (3.48) 

A classical solution u G Xq of the problem (V) is also a weak solution. 

Theorem 3.3. Assume that 

i) the constitutive coefficients satisfy /r > 0, 3 A + 2/r > 0, oi > 0, 02 > 0; 

ii) the loads satisfy the regularity condition f G 

Then there exists one and only one solution of the problem p.48l) . 

Proof. Let us first consider the case 0:2 > 0. The Cauchy-Schwarz inequality, the inequalities (a ± b)^ < 
2 (a^ + b^) and the assumption upon the constitutive coefficients lead to 


(m, v) < C 


II sym Vup -I- II Curl (sym Vu)|p I dv 


II symV'clP -I- II Curl(symVn)|p Idu 


(3.49) 


< C'lklla’o \\w\\xo, 


which means that (•, •) is bounded. On the other hand, we have 


(m,u) = / 

Jn 


2p. II sym Vm|P -I- A [tr(VM)]^ 

-I- pL\ [2ai II devsymCurl(sym Vu)|p -b 2a2 || skew Curl (symVit) 


dv, 


for all u G Xq. Moreover, as a consequence of the properties i) of the constitutive coefficients we have that there 
exists the positive constant c 

{u,u)>c J ^11 sym Vu|p-b II Curl (sym Vu)|p^ dz;. (3.50) 

From linearized elasticity we have the first Korn’s inequality m, that is 

l|VM||L 2 (n) < CII sym VM||z, 2 (n), (3.51) 

for all functions u G 1 / 0 ( 0 ; F) with some constant C > 0, for bounding the deformation of an elastic medium 
in terms of the symmetric strains. Hence, using the Korn’s inequality (13.511) . it results that there is a positive 
constant C such that 

{u,u) > c ^^||Vuf -b II Cur^symyzi)!!^^ dv = c\\u\\l.^. (3.52) 

Therefore our bilinear form (•, •) is coercive. The Cauchy-Schwarz inequality and the Poincare-inequality imply 
that the linear operator If) is bounded. By the Lax-Milgram theorem it follows that (13.481) has one and only 
one solution. The proof is complete in the case 02 > 0. 

Now, we consider the case 0:2 = 0. Using Lemma 15^ it follows that the bilinear form (•,•) is also coercive 
for (22 = 0. Using similar estimates as above the existence follows also in this case and the proof is complete. □ 

Remark 3.4. The Lax-Milgram theorem used in the proof of the previous theorem also offers a continuous 
dependence result on the load f. Moreover, the weak solution u minimizes on Xq the energy functional 


I{u) 


2^11 sym Vzi|p -b A[tr(VM)]^ 


p-L^flaiW devsym Curl (sym Vu)|p -b 2 a 2 || skew Curl (sym Vw)!!^] — (/, u) 


dv. 
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Let us consider v G C“(ri;r) and u a solution of problem (I3.36I1 - (I3.38I1 . Then we obtain 


{u,v) = / {f,v)dv+ / ((cr + symCurl (to)), n — V[(symM) T] : T, u)(ia 

Jn Jon 


(3.53) 


+ / ((1 — n 0 n)(symM).n, (1 — n C) n)(Vr’).n)c?a + / (([symM]'*" — [symM] ).i>,v)ds 


/an 


/ar 


= / / ((cr + sym Curl (to)), n, a)fia — / ^ (to^ x n, {symVv)i)da. 

Jn Jan Jan 

Therefore, the corresponding existence results assures that there exists the weak solution u minimizing on 
C“(n;r) the energy functional 

I[u) = [ 

Jn . 


2^11 sym Vm||^ + A[tr(VM)]^ 

+ [2ai II dev sym Curl (symVM)|p + 2a2 || skew Curl (symVu)||^ 


~ / ((o'+ sym Curl (m)). n — V[(symM) T] : r, M)(ia 

Jan 


dv — (/, u)dv 

Jn 


(3.54) 


/an 


((1 — n 0 n)(symM).n, (1 — n 0 n)(Vw).n)da — / (([symM]'*" — [symM] ).v,u)ds 


/ar 


/n L 


2 /X II sym Vm|P + A[tr(Vu)]^ 


+ /xLg [2 q;i II dev sym Curl (symVM)|p + 2a2 || skew Curl (symVu)||^ 

~ / ((o'+ sym Curl (to)), n, w)(ia + / y {fhj x n, {sym\7v)i)da. 

Jan Jan 


dv — (/, u)(ia 

Jn 


3.3 Traction boundary condition in the Curl (sym Vw)-formulation 
versus the V[axl(skewVM)]-formulation 

In this section we compare the possible traction boundary conditions in the V[axl(skewVM)]-formulation and 
the Curl (sym Vit)-formulation. The conclusion is summarized in Figure Hand Figure H Prescribing Su and 
(1 —n0n).curlu on the boundary means that we have prescribed independent geometrical boundary conditions, 
this is also the argumentation of Mindlin and Tiersten [72], Koiter m, Sokolowski nnni, etc. However, the 
prescribed traction conditions remain not independent, in the sense that g leads to a further energetic conjugate, 
besides t, of u. From this reason we claim that, in order to prescribe independent geometric boundary conditions 
and their corresponding completely independent energetic conjugate (traction boundary conditions), we have 
to prescribe u and (1 — n 0 n)Vu.n. In other words, we prescribe 

/ {t,u)da+ / (g, {1 — n (Si n)Vu.n)da, (3.55) 

Jan Jan 

in which now u and (1 —n0n)Vu.n are independent and g does not produce work against u, see |6T] for further 
detailed explanations. This type of independent boundary conditions are also correctly considered already by 
Bleustein Ml but for the full strain gradient elasticity case only. In order to have a complete overview on the 
subject, in Table I we also summarize the equivalent form of the equilibrium equations. 
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' 6u 

and 

VSu X n 

cannot independently be prescribed 


6u 

and 

V(5u. T 

cannot independently be prescribed 


Su 

and 

(1 — n C) n) S7Su 

cannot independently be prescribed 


Su 

and 

VSu. n 

can be independently prescribed (3 be) 


6u 

and 

(1 — n (g) n) VSu. n 

can be independently prescribed (2 be) 


6u 

and 

(curlJu, n) 

cannot independently be prescribed 

J 

' Su 

and 

(curl Su, t) 

can be independently prescribed (2 be) 

u 

_ Su 

and 

(1 — n 0 n).curl Su 

can be independently prescribed (2 bc)y' 


Figure 3: The possible independent geometrically boundary conditions. The joining bracket means that the 
conditions are equivalent. 


Table 1. Euler-Lagrange equations in various formulations 


Euler-Lagrange equations 
in direct tensor format 

Euler-Lagrange equations 
in indices 

Euler-Lagrange equations for Curl (sym Vtt) 

Div((7 -1- t) -1- / = 0 

tr = Dsym VuVflin(sym Vu) e Sym(3) 

r = sym(Curl fh) S Sym(3) 

m = L>Curl (sym VulfTcurv(Curl (sym Vii)), second order 

Euler-Lagrange equations for Curl (sym Vu) 

{aij + fi = 0 

Oij = jWiin(|(uiq + G Sym(3) 

W = \ + ^jlkffiik,i) e Sym(3) 

Euler-Lagrange equations for V[axlskew Vu] 

Div((j — r)-\-f = 0 

a = DsymVuWiini^ymVu) G Sym(3) 

r = anti(Div m) G so(3) 

m = -Dv[axl(skew Vtt)]^curv(V[axl(skew Vw)]), second order 

Euler-Lagrange equations for V[axl skew Vu] 

{pij ~ H" /i = 0 

.)Wiin(|(«i,f +“j,i)) e Sym(3) 

'^ij — '2^jik'^kl,l ^ 50 ( 3 ) 

Euler-Lagrange equations for V[axlskew Vu] (3rd order) 

Div(f7 — r)-\-f = 0 

a = DsymVuWiini^yrnVu) G Sym(3) 
r = Divm = ^ anti(Divrn) G so(3) 

m = DvVu[^curv(V[axl(skew Vii)])], third order 

Euler-Lagrange equations for V[axlskew Vu] (3rd order) 

{^ij ~ ^ijk,k) J /i “ 0 

aij = -Di Wiin( 2 (^dj ^ Sym(3) 

Tij = “ 2^jik'f^kl,l ^ S0(3) 

mijk = Du ijf,WcuTv ((epfeUfej) „,) 


We outline that there exists a relation between the allowed traction boundary conditions in the Curl (sym Vu)- 
formulation and those from the V[axl(skewVu)]-formulation which we take from |61j 
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~ / {{<7 — T).n,Su) da — / {m.n,ax\{skewVSu)) da 

JdQ, J dQ 


(3.56) 


axl-formulation 


= ~ / {{cr+ T).n,Su)da— / ^ {rhj x n, (sym da . 

Jdn Jon 


Curl-formulation 

or after splitting up with use of the surface divergence theorem on both sides 


/an'- 


{{a + symCurl {m)).n — V[(symM) (1 — n 0 n)] : (1 — n 0 n), 5u)da 


(3.57) 


( 1 ) 


(a) 

— / {{1 — n®n){syTa.M).n,V5u.n)da— / (([symM]'*" — [symM]“). du)ds 


/an 


( 2 ) 


/ar 


(3) 


(&) 


(c) 


/■ 1 1 

= — / ((cr — - antiDiv[TO]). n —-V[(anti(m. n)) (1 — n (8) n)] : (1 — n (8) n), du) da 
dan 2 _2__^ 


(!') 


(o') 


— / (-(1 — n (8) n) anti(m. n).n, Vdu.n)da — f (-([anti(m. n)]^ — [anti(TO. n)] ).i',5u)ds, 
dan 2__^ dar 2__^ 


( 2 ') 


(3') 


ib'] 


(c') 


for all variations Su G C°°(n). Naively, we might expect that the quantities involved have to be equal term 
by term, i.e. (1) = (!'), (2) = (2'), (3) = (3') or (a) = (o'), (b) = (6'), (c) = (c'). However, this is not true, see 
Appendix 

3.4 Principle of virtual work in the indeterminate couple stress model 

3.4.1 Principle of virtual work in Cosserat theory 

Let us first recall that in the Cosserat theory with independent fields of displacement and microrotations the 
internal energy has the form W (Vu, A, VA) in which u : fi —>■ is the displacement and A : H —>■ so(3) is the 
inhnitesimal microrotation. The virtual work principle of the Cosserat theory is given by 


^int _ ^ext 


(3.58) 


where 


•^int _ 


dt 


lT((Vu + tVSu), A + tSA, V(A + tSA)) 


t=o 


P®’'* = J { f ,u)dv + J { axl(9Jl) ,axl(A))da 

body force body couple 

+ / ( t ,u)da+ { axl(0) , axl(A)) da, 

Jan\T , dan\r '—v—' 

surface tractions surface couple 


(3.59) 
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with / : srj \ r C R3 t : af} \ r C R3 -)> ]R3, : rj C IR3 ^ so(3), © : 17 C R3 so(3). From this virtual 

work principle, one obtains the equilibrium equations 


Divcr + / = 0, (3.60) 

Divm + axl(skew cr) + axl(97t) = 0, 

where cr = I7symVukF(sym Vu, A) and m = D^-jWiyu, A,yA), and the boundary conditions 

a.n = t, m. n = axl(0). (3.61) 

In order to obtain these equilibrium equations and the form of the boundary conditions, we have used the fact 
that u and A are independent constitutive variables. 


3.4.2 Virtual work principle in the Curl (sym(Vu))-formulation of the indeterminate couple stress 
model 

Using again the fact that u and 2 curlu = axl(skewVu) ~ A are not independent constitutive variables, and the 
identity 

V[axl(skew Vu)] = [Curl (sym Vu)]^ 

we consider the energy lU(sym(Vu), Curl (sym(Vu))) and the following new form of the virtual work principle 

'pint _ peyit 


(3.62) 


where 


1 p 

P'"* = — / lU(sym(Vit + fVJu), Curl (sym(Vu + tVJu))) 
/ ( / ,u)dv + / ( M ,symVu)du 

Jq. Jq , , 


t=o 


pext _ 


(3.63) 


body forces 


body couple 


/ _( t ,u)da+ / _( 

Jdn\r , Jdn\r 

' surface traction ' 


9 


surface double tractions 

K^, : 17 C 


, {1 — n (Si n)\7u.n)da + (p ( ,v)ds, 

^ _ 1_1 ._ r _ 


with f : n <Z W ^ t : dn\T <Z W ^ 

TT : dr C ]^3_ 

From this virtual work principle, we obtain the equilibrium equations 


edge line force 

so(3), g : 517 \ F C K.^ —>• and 


Div (cr + T + Al) + Div A1 + / = 0, 


(3.64) 


the total body force 


where 


cr = I7symVubF(sym Vu, Curl (sym Vu)), 

T = sym[Curl m] S Sym(3), 
fh = I7curi(symV«)W^(syniVu, Curl (symVu)) 

and the following traction boundary conditions 
[{a + T).n — V[(symM) (1 — n 0 n)] : (1 — n 0 n)] (x) = t{x) — M.n, 

the total traction condition 

[(1 — n 0 n)(symM).n] (x) = [(1 — n 0 n) g] (x), 

{([symM]+— [symM]“). i/} (x) = tt, 


> X G 517 \ F 


X G 5r 


(3.65) 


(3 be) 
(2 be) 

(3.66) 
(3 be) 


where F is an arbitrary open subset of 517. 

In Subsection 13.31 and Appendix El we show that the traction boundary conditions in the Curl (sym(Vit))- 
formulation do not coincide pointwise with those arising from V[axl(skew VM)]-formulation, even if the boundary 
virtual power works are identical. 
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Incomplete-standard boundary conditions in the V[axl(skew Vu)]— formulation |72| 

Geometric (essential) boundary conditions (3-1-2) [correct] 

(1 — n(8) n).curlii|p = (1 — n (gi n).curlu^ £ 

Mechanical (traction) boundary conditions (3-^2) [erroneous] 

((cr - T).n - in X V[(n, (symrn).n>]) = i", 

(1 — n (g) n) m. = (1 — n (g n) /i 

Boundary virtual work 

— fQ^{{(7 — t). n, Su) da — (m. n, axl(skew V5u)) da = 0 ^ 

— (|((T — T).n — X V[(n, (sym m).n)] |, Su) da — {rh.n, |(1 — n g n) [axl(skew V5ti)] |) da = 0 


, or (1 — n g n)Va.n|p = (1 — n g n) Vu^.n £ 
T = Divm = ^ anti(Divm) £ so(3) 


3 be 
2 be 


normal curvature 


$ 


Incomplete-standard boundary conditions in the V[axl(skew Va^]—formulation, written in indices [58| 
Geometric (essential) boundary conditions (3-1-2) [correct] 

or {ui^k^k - Uj^knknjUi) = ^i,k'^k - u^^f^nkrijUi 
Mechanical (traction) boundary conditions (3-^2) [erroneous] 


dn\r ~ 


{^‘^ik'^k ‘^jk^k^j^i^ l0r2\r — hjTljTLi 

Boundary virtual work 

— ii^ij — ^ij) '^j) axl(skew V(5a)) da = 0 


'^ij — G S0(3) 


3 be 
2 be 




fdCl 2 XO 2 fdCt (^^ikl^'^l,k ^jkl^^l,k^j^i) — d 


Figure 4: The incompletely boundary conditions in the V[axl(skew Vu)]-formulation which have been employed 
hitherto by all authors to our knowledge. The virtual displacement is denoted by Su £ The 

number of traction boundary conditions is correct, but the split into independent variations at the boundary is 
incomplete, as shown in [61] . 
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Boundary conditions in the V[axl(skew Vu)]—formulation in terms of gradient elasticity n 

and third order moment tensors, see [61] 

Geometric (essential) boundary conditions (3+2) 

(1 — n[8) n)Vu.n|p = (1 — n Gi n) E or (1 — n (g) n).curlu|p = (1 — n gi n). curin'^ E 

Mechanical (traction) boundary conditions (3+2) 

((cr — Div m). n — V[(m. n) (1 — n g n)] : (1 — n g n)) = t, m = DvV'u[^curv(V[axl(skew Vu)])] 3 be 

(1 — n g n)[m. — (1 ~ “S) ^^) ^ 2 be 

([m. n]"*" — [m. n]“). ^ “edge line force” on dP 3 be 

Boundary virtual work 

— — Divm). n, 5ii) da — fQ^{xn.n,S/Su) da = 0 

— f^Q((a- — Divm). n — V[(m. n) (1 — n g n)] : (1 — n g n), du) da — — n g n)[m. n].n, VSu.n)da 

— /Qp(([Tn- n]^ — [m. n]~).iy, 5u)ds = 0 




equivalent 


Boundary conditions (3+2) in the V[axl(skew Vu)]— formulation in terms of gradient elasticity, third order 
moment tensors, and written in indices, see EH 

Geometric (essential) boundary conditions (3+2) 

"Uj I p — ^j,k^j ^k^i ) | p — ’ 

i^ikl'^l,k '^'i) I p “ 

Mechanical (traction) boundary conditions (3+2) 

^ijk,k') '^j (j^ipk '^k ^ijk '^k'^j'^p') ^ {^^ph I^Q^p ^ijk B+urv (^{^^ijk^k,j') 

(iTlijp Tlj ^pjk '^k^j'^i') ~ hpnpTii 

{[mpjk n^]'^ - [mpjk nk]~). Uj |^p = ^p “edge line force” on dV 


3 be 

2 be 

3 be 


Figure 5: The complete-standard boundary conditions in the V[axl(skew V?x)]-formulation in terms of a third 
order couple stress tensor coming from full gradient elasticity. The virtual displacement is denoted by 5u E 
The summation convention was used in index notations. 
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4 Relation to the Cosserat-micropolar and micromorphic model 


We have seen that it is irrelevant whether we take V[axl(skew Vm)] or Curl (sym Vu) as basic curvature measures 
for the indeterminate couple stress model as long as consistent requirements on F = dQ are considered and the 
following Dirichlet conditions are used both together 

m|p = uo, (1 — n 0 n)Vti.n|p = (1 — n (g) n)VMo-^ ^ '*^Ip = (1 — n 0 n).curlu|p = (1 — n 0 n).curluo- 

The difference of the formulation appears only when considering mixed Dirichlet-Neumann boundary conditions. 
However, when we want to switch from a 4th-order (gradient elastic) problem to a 2nd-order micromorphic model 
or Cosserat model [m EH EH], we need to introduce new independent variables and decide about the useful 
coupling conditions in terms of adding a penalty term. It is also clear that adding more variables it depends on 
the number of the added helds whether the new formulation is weaker softer in the language of a finite element 
context. In general, more degrees of freedom mean weaker response, at the prize of needing to specify more 
boundary conditions. 

We discuss the following cases: 

i) [Cosserat] skewVtt A € so(3). In the case V[axl(skew Vu)] we are led to introduce a skew-symmetric 
variable A G so(3) instead of skew Vu, thus using the curvature tensor V axl(H) together with the coupling 

II skew(VM) — A|p = ^ Ijcurlu — 2axl(A)|p, 

leading to the classical Cosserat model, with a new penalty parameter Hc > 0 known as the Cosserat 
couple modulus. To be more precise, the corresponding minimization problem becomes 


I{u,A) = 


^ ||symVu|p -I- — [tr(symVM)]^ -f Hc\\ skew(Vu) — A\\‘^ 

+ (oill devsym V axl(H)|p -|- tr[V axl(H)]^ -|- a2||skewV axl(H)|p) 


dV 


w.r.t u G iFg(H), A G IFq)!!). In this case, the force-stress tensor is clearly non-symmetric 

a = 2 ^ sym Vu + 2fj,c (skewVu — H) -|- A tr(Vu) 1 ^ Sym(3), (4.1) 

and the couple stress tensor (hyperstress tensor) is given by 

fh = 2fj,L^ (oi dev sym V axl(H) + ai tr[V axl(H)] 1-1-02 skewV axl(H)), (4.2) 


which is also in general non-symmetric. Note that m has now 3 independent length scale parameters. 

ii) [microstrain] sym Vtt i—)■ £ 6 Sym(3). In the case of starting with the representation Curl (sym Vu) 
we are led to introduce a symmetric tensor variable e G Sym(3) instead of sym Vu, thus using the curvature 
measure Curl e together with the coupling 


leading to a “microstrain” theory 


/(u,e) = f 
Jn 


xr^ll sym Vu — e||^, 

, the minimization problem is now 
A 


fj, ||symVu||^ -|- ^ [tr(symVu)]^ -I- sym Vu — e||^ 

devsym Curl ejp-|-/SaljskewCurl sjp) dV —?> min. 


w.r.t u G iJg (H), e G iJg (Curl; H), and, in this case, the force-stress tensor is symmetric 
a = 2fi sym Vu -f 2 (sym Vu — e) -I- A tr(Vu) 1 G Sym(3), 


(4.3) 
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Correct boundary conditions in the V[axl(skewVtt)]-formulation (V(curlti)-formulation), see |61} ^ 

Geometric (essential) boundary conditions (3+2) 

= (1 — n[8) n)V'U.n|p = (1 — n (gi n) .n E or (1 — n [g n).curl'u|j, = (1 — n g n). curin'^ E 

Mechanical (traction) boundary conditions (3+2) 

((fj — t). n — X V[{n, (symm).n)] — ■iV[(anti[(l — n g n)rh. n])(l — n g n)] : (1 — n g n)) 3 be 

(1 — n g n) anti[(l — n g n)rh. = (1 — n g n) h 2 be 

“edge line force” on dF: {([antilm. n]]"*" — [anti[m. n]]“). = tt 3 be 

Boundary virtual work 

— fQ^{{cr — t). n, Su) da — (fh. n, axl(skew V<5ti)} do. = 0 

. n — X V[(n, (sym m).n)] j-, du) da + Jq^ {m.n, | (1 — n g n) [axl(skew V5ti)] |) da = 0 

normal curvature 

— /qq{(o' — r). n — X V[{n, (symfri).n}] — ^ V[(anti[(l — n g n)rh. n]) (1 — n g n)] : (1 — n g n), 6u) da 

~h — n g n) anti[(l — n g n)m. n].n, ^Su.n)da — ^ jQp(([anti[rh. n]]"*" — [anti[m. n]]~). i/, Su)ds = 0 J 


$ 


equivalent 


Alternative equivalent correct boundary conditions in the V[axl(skewVu)]-formulation, see m 

Geometric (essential) boundary conditions (3+2) 

(1 — n g n)Va.n|p = (1 — n g n) Vu^.n E or (1 — n g n).curlii|j, = (1 — n g n).curla^ E 

Mechanical (traction) boundary conditions (3+2) 

((a — r). n — ^V[(anti(m. n)) (1 — n g n)] : (1 — n g n)) j^^N^p = i, 

(1 — n g n) anti[rh. n] = (1 — n g n) /i 

([anti[m. n]]"*" — [anti[m. n]]“). ^^j^p = ^ “edge line force” on dF 

Boundary virtual work 

— fQ^{((7 — T).n,Su)da — (m. n, axl(skew V5a)) da = 0 ^ 

— /^^((o' — r). n — i V[(anti(m. n)) (1 — n g n)] : (1 — n g n), 5u) da 

~h — n g n) anti(m. n).n, 'V5u.n)da — ^ J'gp(([anti[m. n.]]+ — [anti[m. n]]”). z/, 6u)ds = 0 


3 be 

2 be 

3 be 




equivalent 


Alternative equivalent correct boundary conditions in the V[axl(skewVu)]-formulation, written 

Geometric (essential) boundary conditions (3+2) 


in indict 




or {ui^k^k - 'tJ‘j,knk'f^jrii) Ip = ^i^k'^k - u^^^UkTijUi 
Mechanical (traction) boundary conditions (3+2) 

2i^ihk'^ks'f^s ^ijk'f^ks’^s'ajTlfi) ,p{Shp f^h'^p)') |^Q\p “ 

{^ipk'^ks'^s hpnpTli, 

([^ipk^ksns]'^ - [eipk'mksns]~) i^pj^p = “edge line force” on dF 


3 be 

2 be 

3 be 


Figure 6: The possible boundary conditions in the V[axl(skew Vu)] and Curl (sym Vri)-formulation. The equiv¬ 
alence of the geometric boundary condition is clear. The virtual displacement is denoted by Su G 
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m = 7712 


Correct boundary conditions in the Curl (sym Vu)—formulation 

Geometric (essential) boundary conditions (3+2) 

■u|^ = ^ ^ (1 — n (g) n)Vn.n|p = (1 — n g) n) E or (1 — n g) n).curln|j, = (1 — n g) n) curin'^ E ] 

Mechanical (traction) boundary conditions (3+2) / rh x n \ / rh \ 

(cr + r). n — V[(sym M) (1 — n g) n)] : (1 — n g) = t, M = j rh 2 X n _ 1 m = j rh 2 _j 3 

(1 — n g n)(symM).n|^^^P = (1 — n g n) h, ^ ^ ^ ^ ^ ^ 2 

([symM]"*" — [symAf]“). “edge line force” on dF 3 

Boundary virtual work 

^ 3 _ 

— {(a + r). n, Su)da — {rhi x n, (sym V(5n)i) da = 0 

i=l 

— Jq^ {{a + r). n — V[(sym M) (1 — n g n)] : (1 — n g n), <5ii}da — — n g n)(sym M).n, {\/Su).n)da 

— /^p(([sym M]"*" — [symM]“). Su)ds = 0 




Correct boundary conditions in the Curl (sym Vu)—formulation, written in indices ' 

Geometric (essential) boundary conditions (3+2) 

■^■ilp ~ ^ ^jkl‘^l,k'^j'^i') |p “ 

or - Uj^k^k^jUi) Ip = ^^k^k - nP 

Mechanical (traction) boundary conditions (3+2) 

(cTjj + Tij ) 71 j i^ikl'kf^hk'^l + ^hkl'^ik'^l ^jkl'^f^ik'^l'^j'^h ^ikl'^^jk'^l'^j'^h) ,p i^hp '^h'^p) \ qq\^Y “ ^ 

2 {.^ikl'kapkT^l + EpklTTlik'ki'l ^jkl'kO-pkTll 71 j Tli €pkl771jk71l71j71i ) Tip | ^Q'^p — hpTlpTli , 2 bc 

^([^pkimikTii + eikimpkni]'^ - [epkirdikTii + eiki7hpk7ii]~) i^p\q^ = 3 bc , 


Figure 7: The possible boundary conditions in the V[axl(skew Vu)] and Curl (sym Vu)-formulation. The equiv¬ 
alence of the geometric boundary condition is clear. The virtual displacement is denoted by Su G IR.^). 
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and the hyperstress-tensor is given by 

TO = 2 /r (/3i dev sym Curl e + P'i skew Curl e), (4.4) 

which is non-symmetric in general, depending on the material parameters. Note again that tr(Curl e) = 0, 
thus the spherical part of the hyperstress tensor vanishes and to features only 2 independent length scale 
parameters. 

hi) [micromorphic] Vm p. In this case we may introduce a tensor p G instead of Vu, and use the 

coupling 

— p\\^ 

leading to a micromorphic theory |281181j , the minimization problem being 


/(u,p) = f 
Jn 


p ||symVM|p -|- — [tr(symVu)]^ -I- :n'^\\Vu — p||^ 

-I- /r ( 71 11 dev sym Curl (sym p) | P -I- 73 11 skew Curl (sym p) 11 ^) 


dV —^ min. 


w.r.t u G i?g(n), p G i/Q(Curl; 17). We also point out that the force-stress tensor in this formulation will 
be non-symmetric 

a = 2 p sym Vm -I- 2 x’*' (Vu — P) -I- A tr(Vu) 1 (4.5) 

= 2p sym Vm -\-2>c^ skew( Vu — P) + 2 sym( Vu — P) -I- A tr(Vu) 1 ^ Sym(3) 


and the hyperstress tensor (non-symmetric) is given by 

fh = 2 pL\ (71 dev sym Curl (symp) -|- 73 skew Curl (symp)). (4.6) 

In this formulation tr(Curl symp) does not appear since tr(Curl symp) = 0, Vp S R^^^. Thus the spherical 
part of the hyperstress tensor vanishes, and to features only 2 independent length scale parameters. 


iv) [relaLxed micromorphic] for comparison with other extended continuum models we present the rel^Lxed 
micromorphic model [81] [35] |63l [62] [79] . In the relaxed micromorphic model, the minimization problem 
is of the type 


I{u,p) = [ 
Jn 


p ||symVM||^ -I- — [tr(symVu)]^ -I- || sym(VM — p)|p 

-I- L^{oii II dev sym Curl p|p -|- 02 ||skewCurl p|p -|- 03 [tr(Curl p)]^ 


dV 


w.r.t u G IlQ(fl), p G Pq (Curl; 17), and the corresponding force-stress tensor is symmetric 

a = 2p sym Vm + 2 sym(Vu — p) -I- A tr(VM) 1 G Sym(3) (4.7) 

and the hyperstress tensor is given by 

fh = 2pLl (oi dev sym Curl p -I- 02 skew Curl p -|- 03 tr[Curl p] l), (4.8) 


with a non-vanishing spherical part of the hyperstress tensor. Note that to has 3 independent material 
parameters. 


v) we have also proposed a further relaxed micromorphic model [81] [35] (63] [62] [^ , in which case the 
minimization problem is of the type 


/(m,p) = [ 
Jn 


p ||symVu|p -I- ^ [tr(symVM)]^ -I- xr”*" || sym(Vu — p)|p 


-I- pL^[ai II dev sym Curl p||^ -|- ||skewCurl p|p) 


dV 
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(4.9) 


w.r.t u € 77g(ri), p G i/g (Curl; fJ), the corresponding force-stress tensor is symmetric 

a = 2 p sym Vu ■^2yi^ sym(VM — p) -I- A tr(Vu) 1 G Sym(3), 

and the hyperstress fh is trace free 

m = 2pL^(ai dev Curl p-f 02 skew Curl p). (4-10) 

The further relaxed micromorphic model remains well-posed |35) . A still weaker variant is v) with a 2 = 0. 
Whether this choice is mathematically well-posed is yet unclear. 


5 Conclusion 

Our new symmetric-conformal Curl (sym Vw)-reformulation has the following crucial properties setting it apart 
from existing formulations of couple-stress models: 

• the local and the nonlocal force stress tensors (cr, r) are both symmetric, while the couple stress tensor 
m is symmetric in the conformally-invariant model. 

• the curvature energy is conformally invariant and the couple stress tensor fh vanishes for conformal 
displacement. 

• the model has only one additional length scale parameter, similar to the modified couple stress model. 

• the model is derived with consistent boundary conditions: either 5 geometrical conditions or 5 mechanical 
(traction) conditions. The mechanical conditions are separated into force stress tractions and couple stress 
tractions and correspond to completely independent boundary conditions. 

• for mixed Dirichlet-Neumann boundary conditions the model does not reduce to the modified indetermi¬ 
nate couple stress model. 

The energies in both possible formulations (in terms of V[axl(skew Vrt)] or Curl(symVu)) are the same, 
differences appear only once traction boundary conditions are specified. The need for prescribing this or that 
boundary conditions determines which model should be used. 

In a polar gradient elasticity model we could influence directly continuum rotations without prescribing wlp = 
0. But this should only be possible in a theory which extends beyond mechanics: for example to magnetic or 
electric effects, i.e. needed for particular loading and boundary conditions which excite particular micro-rotations 
(“polarization”). In contrast, in a non-polar elasticity model it is not possible to influence directly continuum 
rotations but a non-polar model is applicable and much more appropriate in a purely mechanical context (see 
Figure Hi]). The case iv) in Fig. ITOl needs mathematical discussion. The extension of the well-posedness to the 
finite strain case in which the corresponding Lagrangian may be written as IF = W{U) + IFcurv(b^, Curl U), 
where F = RU is the polar decomposition is yet missing. Some steps in this direction are presented in |54j . 

6 Epilogue: Much ado about nothing 

We have seen how much effort it took us to derive the consistent boundary conditions in the indeterminate couple 
stress model. The conceptual advantage of not having to discuss the physical meaning of independent degrees 
of freedom is, now, more than outweighed by the burdensome interpretation of traction boundary conditions. 
Nevertheless, all presented formulations are shown to be mathematically well-posed. In the last part of the 
paper we have had a look at 2nd-order (micromorphic) approximations of the given gradient elastic models. 
In these micromorphic models, the boundary conditions are completely transparent. However, it seems that in 
this larger class of models there is yet another variant (the relaxed micromorphic model with integral boundary 
coupling) which combines conceptual simplicity, symmetry of force stress tensor and symmetry of moment stress 
tensor, simplicity of traction boundary conditions and well-posedness to make it superior to all other presented 
formulations. With hindsight, we understand why the indeterminate couple stress model had been abandoned 
in the late ’60ies. For us it is a mystery how it was possible at all to identify material parameters in a theory 
in which boundary conditions had not been conclusively settled? 
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3 -f- 3 = 6 geometric be: u\ 
3 -f- 3 = 6 traction be: cr.nl 


r’ 

an\r’ 


d(A)| 


r 

rh.nl 


dn\r 


ii) sym Vu Sym(3), microstrain model 

G 2 ^ sym V-u + 2 sym(Vii — e) € Sym(3) 
rh = 2 Curl e E 0t(3) 

8 ^ ||symVii||^ + 1| sym Vu — ^P + ^L^|| Curl 

invariant under: u u + W.x + 6, 6 E 

s I—^ W E so(3) 

well-posed: -u E e E//(Curl; 12) 


3 -f- 3 • 2 = 9 geometric be: u\ 
3 -f- 3 • 2 = 9 traction be: a.n\ 


ir’ 

dn\T^ 


^ir 
m. Ta I 


dn\r 


/^ii) Vn. I 


p ^ Sym(3), micromorphic > 

a 2 p sym Vu + 2 x”^ (Vu — p) ^ Sym(3) 
m = 2 /i Curl symp E 0l(3) 

8 ^ p ||symVu|p + x+|| Vu — p|P + /i || Curl symp|p 

invariant under: u + W.x + b, 6 E 

p !->■ p + VC, VC E so(3) 

well-posed: u E p E -L^(S2), symp E //(Curl; 12) 

3 -h 3 • 2 = 9 geometric be: u|p, [symp].TQ;|p 


skewVii > A E so (3), degenerate Cosserat model 
a 2 p sym Vu -|- 2 skew(V'U — A) ^ Sym(3) 
m = 2 p,L^ skewV axl(A) E so (3) 

8 p, ||symVu|p -h pc\\ skew Vu — A|p H- pL^ ||skewVaxl(A)|p 

invariant under: u u + W.x + b, b £ 

Ah^A-vTC, TyEso(3) 
well-posed: u E //^(12), axl(A) E //(curl; 12), 
degenerate curvature energy, 

3 + 2 = 5 geometric be: u|p, (axl(A), Tci}|p 

2 = 5 traction be: (m.n, Tq) [^^^P_ ^ 


i) skewVu e->- A E so(3), Cosserat model 
G 2 p sym Vu -\-2 pc skew(Vu — A) 0 Sym(3) 
in = 2 pL^\/ axl(A) E 0l(3) 

8 p ||symVu|p + pc\\ skew Vu — A|p+ pL^ ||Vaxl(A)|pl 
invariant under: u \-y u + W.x + 6, 6 E R^ 

Ah-^A + ly, TCeso(3) 
well-posed: u E //^(12), axl(A) E //^(12) 



constrain 
Vu = A 

Pc ^ OO 


constrain 
Vu = p. 



Hadjesfandiari-Dargush model ^ 

G 2 p sym Vu E Sym(3) 

in = 2 p skewV[axl(skew Vu.)] E so(3) 

r = 2pL^ anti[Div(skewV[axl(skew Vu)])] E so(3) 

8 p llsymy-ulp pL^ ||skewV[axl(skew Vu)] |p 

invariant under: u u W.x + 6, 6 E R^ 

Vu !->■ Vu + VC, VC E so(3) 
well-posed: u E (12), curl u E //(curl; 12) 

5 geometric be: u|p, (curlu,T q,)| p 
.5 traction be y 


limit model-V axl skew(Vu) formulation ^ 

G 2 p sym Vu E Sym(3) 

in = 2 p V[axl(skew Vu)] = p V(curl u) E 0l(3) 
T = 2pL^ anti[Div (V[axl(skew Vu)])] E so(3) 

8 ^ p ||symVu|p + pL'^ ||V[axl(skew Vu)]|p 

invariant under: u u + W.x + 6, 6 E R^ 

Vu ^ Vu + VC, VC E so(3) 
well-posed: u E //^(12), curlu E //^(12) 

5 geometric be: ulp, (1 — n 0 u) Vu.njp 
5 traction be j 


constrain 
sym Vu = e 


$ 


r = ar2. For Dirichlet + Neumann 

mixed boundary conditions: different models. 


'limit model-Curl (sym Vu) formulation 

(T ~ 2 sym Vn G Sym(3) 

fh = 2 Curl (sym Vn) G 0l(3) 

T = 2 lip sym Curl (Curl (sym Vn)) G Sym(3) 

£ ^ fj, ||symVn|p + /ipW Curl (sym Vn)|p 
invariant under: n i—>■ n + W.x + b, b £ R® 

Vn !->■ Vn + W, W G 3o(3) 
well-posed: n G sym Vn G -ff(Curl;f2) 

5 geometric be: n|p, (1 — n (g) n) Vn.n|p 
V^^raction be___^ 


3 -}- 3 • 2 = 9 traction be: a.n\ 


dn\r’ 


Figure 8: Different possibilities of lifting the variants of the 4th.-order indeterminate couple stress model to 
a 2nd.-order micromorphic or Cosserat-type formulation formulation. In the penalty case, all the considered 
alternatives lead to the same limit model provided only geometric boundary conditions are imposed. It is not 
surprising that the limit model has some peculiarities since the limit procedure is itself singular. Note that 
different micromorphic or Cosserat type formulations generate different sets of boundary conditions. Here Tq, 
q; = 1,2 denote two independent tangential vectors on the boundary. 


’ V[axl(skew Vu)] = [Curl (sym Vm)]^ but 
V axl(skewp) [Curl symp]"^ for p 


Figure 9: Integrability conditions. 
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iv) the relaxed micromorphic model 
cr ~ 2 /i sym Vtx + 2 sym(Vw — p) € Sym(3) 
m = 2 fiL^ Curl p E 0l(3) 

£ ^ p ||symVu|p + || sym(V'U — p)|p + pL^ || Curl p||^ 

invariant under: u-\- W.x + 6, 6 E 

p^p + W, Ty,WEsD{3) 

well-posed: u E p E i^(Curl; n) 

3 -f- 3 • 2 = 9 geometric be: u|p, p. Tq; |p 
3 -|- 3 • 2 = 9 traction be: cy.n \ 


lan\r’ 


constrain 
sym Vti = sym p 


ian\r 


v) the further relaxed micromorphic model ^ 

a ^ 2 p sym Vu -h 2 x"*" sym(Vu — p) E Sym(3) 
m = 2 pL^ dev Curl p E sl(3) 

S p ||symVu|p + x"*" || sym(Vu — p)||^ + p || dev Curl p|p 


limit model: incompatible linear elasticity 

(7 2 p symp E Sym(3) 

m = 2 pL^ Curl p E 0l(3) 

£ ~ M l|symp||2 + ^lLl\\ Curl^ll^ 
invariant under: p p -\- W, W E so (3) 
well-posed: p E -H’(Curl; Q) 

3*2 = 6 geometric be: p. Tq |p 
3*2 = 6 traction be: Curl p. Ta I 


lan\r 


invariant under: u u ^ W.x -|- 6, 6 E B 

p^p-\-W, w^w, 

well-posed: u E p E //(Curl;!!) 


W,W eso{3) 


3 -h 3 * 2 = 9 geometric be: 
3 3 * 2 = 9 traction be: 


Ip, P- 'Tq; Ip 


laQ\r 


vi) another relaxed micromorphic model 

(7 2 p sym Vu -h 2 x"*" sym(Vu. — p) E Sym(3) 

m = 2 p sym Curl p E Sym(3), symmetric 
S p ||symVu|p + x"*" || sym{Vu — p)||^ + p || sym Curl; 
invariant under: u u + W.x + 6, 6 E 

p^p + W, W^w, TF,Weso(3) 

well-posedness not clear : u E p E //(Curl;n) 

3 -f- 3 * 2 = 9 geometric be: u|p, p. Tq |p 
3 -h 3 • 2 = 9 traction be: <j.n j^Q^p, 


/^ii) relaxed micromorphic model with integral boundary coupling 

cr 2 p sym Vu -|- 2 x^ sym(Vii — p) E Sym(3) 
m = 2 pL^ Curl p E 0l(3) 

E p llsyrnyujl^ -h x+ || sym(V'U — p)|p -h pL"^ || Curl p|p 

+ /p ||{Vu — p) X ri|pda 

invariant under: u u W.x -h 6, 6 E 

pn^p + W, W^W, Ty,WES0(3) 

well-posed: u E i/^(!!), p E//(Curl; £!) 

3 geometric be: u|p 



Classical linear elasticity 
<7 2 p sym Vu E Sym(3) 

E p ||sym Vulp 

invariant under: u + W.x + b 

Vu !->■ Vii -I- VF, W E so(3) 
well-posed: u E 'Vu E //(Curl;!!) 

3 geometric be: u|p 
3 traction be: a.n j^Q^p 


Ln —^ OO 


limit model: incompatible linear elastici^\ 
a 2 p symp E Sym(3) 
m = 2 pL^ sym Curl p E 0l(3) 

E p ||symp||^ + pL‘^\\ sym Curl p|p 
invariant under: p p W, W E so (3) 
well-posed: p E //(Curl;!!) 

3 2 = 6 geometric be: p. Tq | p 
3-2 = 6 traction be: m. Ta 


3*2 = 9 traction be: a.n 


dn\r 


TT, m. Ta 


dn\r 


Figure 10: For comparison, the same situation as in Fig. [5] for the relaxed micromorphic model and the further 
relaxed micromorphic model |81j . The first two micromorphic formulations are well-posed. In our view these 
models are advantageous as compared to the models in Fig. [H Note also that the boundary conditions for 
the new microdistortion field p has 6 degrees of freedom, the only physical and mathematical possible choice is 
p X n|p = 0 on F, where Dirichlet-boundary conditions u|p = 0 are prescribed. In the case of non-homogeneous 
(non-zero) boundary prescription u on F, we would need to modify the total energy by adding as a boundary 
term Jp || (Vu — />) x nj^ds. This will introduce a certain coupling at the boundary. The model in this format 
is still well-posed and awaits to be further investigated. A nonlinear modification of this model is investigated 
in [Si]. 
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Not only rotational interaction 
gradient elasticity: Wcutv{D^u) 
m = Dj^ 2 uWcutv{D'^u), third order J 


Size effects 
>’’Smaller is stiffer” 


Strain gradient elasticity 

Wcurv(V (sym Vu)) 


''V (sym Vii) ’ 


/(V (sym Vw)), third orden 


Only rotational interaction through continuum rotation Wcurv(V(curl li)) = Wcurv (Curl (sym Vw)) 


Cauchy-Boltzmann axiom 
^■total = ®’total 
purely mechanical interactit 


k = V[axl(skew Vu)] = -^Vcurlu/ 

^ ^ X 

/ polar gradient elasticity ' 

j indeterminate couple stress 

= D-^Wcurvik), second orde^ 


Grioli-Toupin-Koiter-Mindlin 


fi II sym Vii||^ 

||Vaxl(skew Vu)!!^ 

2 curvature parameters 
couple stress tensor rh G 9l(3) 
nonlocal force-stress r € so (3) 
unique solution 
extendable to a 
Cosserat formulation 
axl(A) G H^{n) 


Modified couple stress modelA 
conformally invariant _ y 

II II sym Vulp 

-h/i II dev sym V[axl(skew Vu)]|p 

1 curvature parameter 

couple stress tensor m E Sym(3) -4 

nonlocal force-stress r E so (3) ^ 

unique solution 

extendable to a 

Cosserat formulation 

axl(A) G H^{n) 


integrability condition 


difference of formulation 
is a Null-Lagrangian in 
the balance of linear momentu 
and sym m = sym fn 


Hadjesfandiari-Dargush 


fi II sym Vu||^ 

II skew V[axl(skew Vu)]|p 
1 curvature parameter 
couple stress tensor rh E so (3) 
nonlocal force-stress r E so (3) 
unique solution 
extendable to a 
Cosserat formulation 
axl(A) E /f(curl; £7) 


: = Curl (sym Vu) 


non-polar gradient elasticity j 
.rh = D^Wcurvik), second orders 


Curl -formulation 


fi II sym Vu|p 

-\-I^L‘^ II Curl (sym Vu)|p 

2 curvature parameters 

couple stress tensor rh E 0l(3) 

nonlocal force-stress r E Sym (3) 

unique solution 

extendable to a 

microstrain formulation 

£ E H{Cur\;n) 


T := anti Divm G so(3) 
r := sym Curl m £ Sym(3) 


: /J {2 dev sym Curl (symVu) 
+ 2 Q 2 skew Curl (sym Vu) } 


: 1 ^ dev sym V [axl(skewVu)] 

+ 02 skew V[axl(skew Vu)]} 


(^new conformally-invariant model^ 

II sym Vtt|p 

+fM II dev sym Curl (sym Vii)|p 
1 curvature parameter 
couple stress tensor m G Sym(3) 
nonlocal force-stress f G Sym(3) 
unique solution, well posed 
extendable to a 
microstrain formulation 
e G H{Cuvl ;Q) 


Figure 11: Some models involving size effects and their interrelation. The given energy expressions are only 
meant to represent the main features of the models: /i||symVu|p represents the linear elastic energy which 
may be anisotropic (C.sym VUjSymVtt), while /iL^||V[axl(skew Vm)]|P represents the curvature energy, which 
could also be anisotropic. 


r = /i sym Curl {2 ai dev sym Curl (sym Vu) + 2 0:2 skew Curl (sym Vm)} = sym Curl (m) € Sym(3) 


T = fiL'l anti Divjai dev sym V axl(skew Vm) -I- 0:2 skew V[axl(skew Vm)]} = antiDiv(m) G so(3) 


Div(r — r) = 0 


Figure 12: The two possibilities of defining a nonlocal force stress tensor: either r is symmetric in the 
Curl (sym Vw) formulation or t is antisymmetric in the V[axl(skew Vm)] formulation. The difference between 
both stresses is a divergence-free stress field (a self-equilibrated force field). 
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A The traction boundary conditions in the Curl (sym Vrt)-formulation 
and in the V[axl(skewVii)]-formulation are different 


In this section we prove the claim from Subsection lO, i.e. we show that the possible traction boundary conditions in the 
V[axl(skewVri)]-formulation and the Curl (sym Vti)-formulation are different. 

We consider a point P at the boundary and we show that (2) ^ (2') in this point. Without confining, we consider that the 
system of coordinates is initially chosen such that the normal vector on the boundary at this point P is n = e\ := (1, 0, 0). Since 
there are no derivatives in ( 2 ) and ( 2 '), it is enough to prove that 


On one hand, we have 
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and therefore 
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On the other hand, we obtain 


anti(m. ei).ei = anti(mii, m 2 i, m 3 i).ei = 
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Let us now point out that m and m are not independent, see Figure fTTI Considering the case Q 2 = 0 we have m = m € Sym(3), 
while considering ai = 0 we have m = —m E so(3). Therefore, in the conformal invariant case a 2 = 0 and also in the case ai = 0, 
since the condition iixn does not hold true, it follows that ( 2 ) 7 ^ ( 2 ^). 

If the boundary conditions imply the continuity of m and m, then (3) = (3^) = 0. However, if m and m are not continuous 
across the curve OF, considering again a point P E ^F and considering, without confining, that the system of coordinates is initially 
chosen such that the normal vector on the boundary at this point P is n = ei := ( 1 , 0 , 0 ) and 1 / = 63 := ( 0 , 0 , 1 ), we prove that 


symM.e 3 7 ^ anti(m. ei).e 3 . 
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We remark that in both particular cases, the conformal invariance model Q 2 = 0 and the case q;i = 0, we deduce that symM.e 3 7 ^ 
anti(m. ei). 63 , in general. However, even if symM.e 3 7 ^ anti(m. ei).e 3 , the jump (c) may coincide with the jump (c'). 

Let us remark that in order to compare (1) and (!') we may not proceed as above. However, we will prove that (1) 7 ^ (H) in a 
specific situation. We assume that there is an open subset to C <9^1 such that on to the normal vector n is constant. Let us consider 


a point P G (jJ. We may assume for simplicity that n = ei at all points 
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Moreover, we obtain 
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Concluding, (1) = (!') if and only if 
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which holds not true, in general. In the conformal case, the above condition reads 
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B From second order couple stress tensors to third order moment 
stress tensors and back 


Let us consider the general anisotropic case and 

Wcurv{D^u) = {C.D^u, D^u)^3x3x3 , Wcurv(Curl (syui Vu)) = (L. Curl (symVu), Curl (sym Vu))]^3x3, 

where 


C = (Cyfc(^„):IR3x3x3^R3x3x3 L = (Li.xi) 

Let us also consider the tensors 

m := Diy 2 ^WcuTv{D‘^u) and m := (sym Vii)M4urv(Curl (sym Vn)), 

which for our anisotropic case are 

m := and m := L. Curl (sym Vu). 


Since 


D^u = Lin(V (sym Vii)) = A.(V (sym Vii)), Uk^ij = £ifc,j - Sjk,i - £ij,k, 
where s = symVu, we obtain 


m := CA.(V (sym Vu)) = B.(V (sym Vu)). 
The next problem is to find particular form of the tensor C, for which we have 

(m, D^77)g3x3x3 = (m, Curl(symV7t)}g3x3, 
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or equivalently 


{CA.(V (sym Vtx)), A.{V (sym Vii)))]g 3 x 3 x 3 = (L- Curl (sym Vu), Curl (sym Vtx))]^ 3 x 3 • 


This is equivalent to 

{A^CA.(V (sym Vu)), (V (sym V'u)))]^3x3x3 = (L. Curl (sym Vu), Curl (sym Vu))]g3x3 . 
Let us denote in the following 


B := A^CA. 


We consider a specific form of the tensor B in terms of another tensor L : —>• R^^^ such that 

{B.(V (sym Vu)), (V (sym Vu)))]g 3 x 3 x 3 = (L- Curl (sym Vu), Curl (sym Vii))]^ 3 x 3 • 

Let us show how to obtain the tensor B if L is given, such that the last identity holds true. We first remark that a tensor 
]g.R3x3x3 ^3x3x3 jg uniquely defined by the fourth order tensors 

Bim = (B>ijkmnp)im, ■ K®''® ->■ (B.l) 

and 
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(B. V(sym Vii), V(sym V'u)}g3x3x3 := (Bi^-V(sym Vu)m, V(sym Vii)i)j{3x3, 

where Einstein’s summation rule is used. Let L : be a given fourth order tensor. We may write this tensor in the 

form 

L = (L\L2,iL®), L‘:R®^R®. 

Let us define the tensor B by EIJ where 

2 skew anti L*’ axl skew for i = m 
0 for i ^ m. 

For the particular form lES of Bjy^, using the formula 2 axl skewV(sym Vu)^ = curl (sym Vu)i, we obtain 

(Bi^^ .V(sym Vu)m, V(sym Vu)i)]g3x3 = 2 (skew anti L’^. axl skew V (sym Vu)^, V(sym \^u)i)^3x3 

= 2 (antiL*. axlskewV(sym Vu)i, skewV(sym Vu)i)]^3x3 
= 4 (L*. axl skewV(sym Vu)^, axlskewV(sym Vu)i)]g 3 
= (L^.curl (sym curl (sym Vu)i)]^ 3 . 

Now, we conclude that for a given fourth order tensor L : R^^^ ^ R^^^ the tensor L defined by i IbU and (HU is such that 

(B. V(symVu), V(symVu))]g3x3x3 = {Bi^.V(sym V(sym Vn)*}j^3x3 (B.4) 

= (L^.curl (sym Vu)i, curl (sym Vii)i)]g3 = (L. Curl (sym Vu), Curl (sym Vu))]^3x3 • 

In conclusion, we have found a tensor C given by 

C := ABA^, 

where B is given by Hdj and such that 

(C.-D^w, D^ii)]^3x3x3 = (L. Curl (sym Vii), Curl (sym Vu))]^3x3, 

or equivalently 

(m, D^n)g3x3x3 = {m,Curl(symVtt))g3x3, 

or equivalently 

WcurviD^u) = Wcurv(Curl (sym Vi/)). 


C The name of the indeterminate couple stress model 

Regarding the name of the indeterminate couple stress model, Paria [911 p. 1] writes: “...it has led to the difficulties that the 
anti-symmetric part of the stress dyadic as well as the isotropic part of the couple-stress dyadic remain indeterminate. These inde- 
terminacies are perhaps due to the fact that the rotation vector, defined above, is not independent but depends on the displacement 
vector”. The theory has a variety of names, such as “Cosserat theory with constrained rotations” (Toupin, 1964), “Couple stress 
theory” (Koiter, 1964), “Indeterminate couple stress theory” (Eringen, 1968), “Cosserat pseudo-continuum” (Nowacki, 1968). Erin- 
gen writes m- “At this time [in the 1960, our addition] also popular was a theory of indeterminate couple stress which is mostly 
abandoned now [1998]. In this theory, the axisymmetric [skew-symmetric] part of the stress tensor is redundant and it remains 
indeterminate”. Schafer m called “indeterminate couple stress model” as pseudo-Cosserat-continuum of the triedre caches (see 
also [89]). 

If the microrotations A G so(3) are constrained to be equal to the macrorotations skewVu, the Cosserat model reduces to the 
couple stress theory. This corresponds to the case //c oo, for which the antisymmetric part of the strain tensor skew(Vu — A) 
and the spherical part of the curvature tensor tr(V axl(skewA)) tend to zero. Consequently, by energetic duality the antisymmetric 
part of the Cauchy stress tensor skew(cr) in the Cosserat model, and the first invariant of the couple stress, namely tr(m) do not 
appear in the formulation of the virtual work principle as well as in the constitutive equations. The first invariant of the couple 
stress remains “indeterminate” and it is taken to be equal to zero m- Now, the skew-symmetric part of the total force stress 
tensor now is not constitutively determined, but can be obtained from balance of momentum. 
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